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1 Introduction 



Nucleoli substructure has been inves- 
tigated through various high-energy experi- 
ments. Electron or muon projectile is ideal 
for probing minute internal structure of the 
nucleon. The reaction is illustrated in Fig. 
|1.1| , where the virtual photon from the lepton 

interacts with the target nucleon. Its cross „ ■ n k . 

y N 

section is related to two structure functions 

Fi and Fa depending on transverse and Ion- Figure 1.1: Lepton-nuclcon scattering. 
gitudinal reactions for the photon. They depend in general on two kinematical vari- 
ables Q^ = —q^ and x = Q'^/2p ■ q where q is the virtual photon momentum and p 
is the nucleon momentum. These structure functions provide important clues to in- 
ternal structure of the nucleon |1|, 0, |^. It is known that the structure functions are 
almost independent of Q^, which is referred to as Bjorken scaling. It indicates that 
the photon scatters on structureless objects, which are called partons. The partons 
are now identified with quarks and gluons. The cross section is calculated by the 
lepton scattering on individual quarks with incoherent impulse approximation, then 
the structure functions are described by quark distributions in the nucleon: for exam- 
ple F2(x, Q^) = J^i^l^ilii^yQ^) + Qii^jQ'^)]- Because the variable x is the light-cone 
momentum fraction carried by the struck quark, the structure function F2 suggests 
quark-momentum distributions in the nucleon. 

Quark-antiquark pairs are created perturbatively according to Quantum Chromo- 
dynamics (QCD) so that there could be infinite number of quarks and antiquarks in 
the nucleon. A meaningful quantity is, for example, the difference between quark and 
antiquark numbers. It is certainly restricted by the baryon number and charge of the 
proton. The valence-quark distribution qy is defined by g^, = q — q, then the quark 
distribution is split into two parts: valence and sea distributions. With the definition 
of the valence quark, the sea-quark distribution is given hj q^ = q — q^ = q. The 
valence quarks are the "net" quarks in the nucleon. On the other hand, the sea quarks 
are thought to be produced mainly in the perturbative process of gluon splitting into 
a qq pair. Because u, d, and s quark masses are fairly small compared with a typical 
energy scale in the deep inelastic scatting, the splitting processes are expected to occur 
almost equally for these quarks. Therefore, it was assumed until rather recently that 
the sea was flavor symmetric {u = d = s) [Q . 

Both the valence and sea contribute to the electron or muon cross section, so that 
other processes have to be used in addition for studying the details of the sea. Valence- 
quark distributions are obtained in neutrino interactions. Once the valence distribu- 
tions are fixed, the antiquark distributions in the nucleon are estimated from electron 
and muon scattering data or independently from Drell-Yan processes. The fiavor- 

1 



symmetric antiquark distributions in u, d, and s had been used for a while; however, 
neutrino induced dimuon events revealed that the strange sea is roughly half of the u- 
quark or d-quark sea P, || . It had been, however, assumed that antiquark distributions 
u and d are same. If they are different, it should appear as a failure of the Gottfried sum 
rule 1^. The sum rule was obtained by integrating the difference between the proton 
and neutron F2 structure functions over x, Iq = J dx{F2 — F2)/x = 1/3. There is an 
important assumption in this sum rule, and it is the light antiquark flavor symmetry 
u = d. If it is not satisfied, the Gottfried sum rule is violated. However, it should be 
noted that the sum rule is not an "exact" one, which can be derived by using current 
algebra without a serious assumption. Therefore, the fundamental theory of strong 
interaction, QCD, is not in danger even if the sum-rule violation is confirmed. 

There is an earlier indication of the sum-rule violation in the data at the Stanford 
Linear Accelerator Center (SLAG) in the 1970's [Jg=0.200± 0.040] [g. The analysis in 
1975 showed a significant deviation from the Gottfried value 1/3. However, no serious 
discussion could be made on the possible violation because the smallest accessible 
X point in the experiment was a;=0.02 and there could be a significant contribution 
to the sum from the smaller x region. Nevertheless, it is interesting to conjecture a 
possible physics mechanism of the sum-rule violation. Because the integral is given by 
Ig = J dx{u + u — d — d) = 1/3 + (2/3) / dx{u — d), the fact that the measured value 
Jg=0.200 is smaller than 1/3 suggests a J excess over u in the nucleon. Proposed ideas 
for creating the flavor asymmetry in the 1970's are, as far as the author is aware, a 
diquark model and a Pauli blocking mechanism |Ty, |n[]. The details of these models 



are discussed in sections [4.2| and [4.4] . Other experimental information came from Drell- 
Yan processes. Fermi National Accelerator Laboratory (Fermilab) E288 Drell-Yan data 
in 1981 [0 suggested also a flavor asymmetric sea: u = d{l—x)^''^^. Later, the sum rule 
was tested by the European Muon GoUaboration (EMG) in 1983 and 1987 [|l^. The 
1987 analysis indicates that the sum is 0.197 ±0.011(stat.) ±0.083{syst.) in the region 
0.02< X <0.8, and the extrapolated value is 0.235]};o;q99. Again, the data suggested a 
significant deficit in the sum rule. It was, however, not strong enough to surprise our 
community because the measured difference was still within the standard deviation. 
The another muon group at the European Organization for Nuclear Research (GERN), 
Bologna-GERN-Dubna-Munich-Saclay (BGDMS) collaboration, also investigated the 
sum rule in muon scattering on the hydrogen and deuterium [JMl . The BGDMS result 



in 1990 is 0.197 ± 0.006(stat.) ± 0m6{syst.) in the region 0.06< x <0.8 at Q^=20 
GeV^. Their estimate of the small x contribution is between 0.07 and 0.22, so that the 
result could be consistent with 1/3. 

Although the sum rule was proposed in 1967, there is little progress in 1970's and 
1980's. The crucial point was, as it is common in most sum rules, the lack of small 
X data with good accuracy. The first clear indication of the sum-rule breaking was 
suggested by the New Muon GoUaboration (NMG) in 1991 []r3|. They obtained data 
with X as small as 0.004 by using a GERN muon beam. They fitted i^|'(x) — F2{x) data 
by a smooth curve and extrapolated it into the unmeasured small x region. According 



to the NMC, the integral Iq became 0.240±0.016, which is approximately 28% smaller 
than the Gottfried sum. Their reanalysis in 1994 indicates a similar value 0.235±0.026. 
Considering the small errors, we conclude that the light antiquark distributions are not 
flavor symmetric and we have a d excess over u in the proton. 

Recent measurements of F2/F2 by the Fermilab-E665 |T^ and the HERMES ||17| 



collaborations agree with the NMC results. Estimate of the Gottfried sum is not 
reported yet; however, the agreement of F2/F2 suggests a violation of the sum. More- 
over, the charged-hadron-production data by the HERMES support the NMC flavor 
asymmetry. 

On the other hand, there are existing Drell-Yan data. As it was mentioned, the 
Fermilab-E288 in 1981 suggested a d excess over u [O]. However, later Fermilab-E772 



collaboration data showed no signiflcant flavor asymmetry [|1^ in 800 GeV proton- 
induced Drell-Yan measurements for the deuteron, carbon, and tungsten. Strictly 
speaking, these nuclear data cannot be compared with the NMC results because no- 
body knows how large nuclear modiflcation is. A possible nuclear modification oi u — d 



is discussed in section 4.6. There are data from p-p and p-d Drell-Yan processes by 



the NA51 collaboration |^ at CERN. The data indicated large flavor asymmetry 
'u/(i=0.51±0.04±0.05 at x=0.18. It is again a clear indication of the flavor asymmetry 
in the light antiquark distributions. In order to get more information for the asym- 
metry, the E866 experiment is in progress at Fermilab by measuring the Drell-Yan 



processes pO[. Preliminary data also indicate u < d which could be consistent with 
the NMC. The existing E288, E772, and NA51 Drell-Yan results and the details of the 
Drell-Yan processes are explained in section ^ together with other processes, which are 
sensitive to the flavor asymmetry u — d. 

Next, we discuss a brief outline of theoretical studies. First, there is a conservative 
view that the Gottfried sum is satisfled without the u/d asymmetry by including a 
signiflcant contribution from the small x region [x < 0.004) |^. However, this idea is 
not consistent with the NA51 data. We also note that perturbative corrections to Iq 
are fairly small and it is of the order of 0.3% at Q'^=A GeV^ |2g, ^, 0, ^. The small 
correction is from the q -^ q splitting process which can occur in the next-to-leading- 
order (NLO) case. If the sum-rule violation or the flavor asymmetry is conflrmed, it 
should be explained by a nonperturbative mechanism. 

A reliable way of treating nonperturbative phenomena is to use lattice QCD. Al- 
though real lattice calculation of the Gottfried sum is not available at this stage, scalar 



matrix elements were evaluated in Ref. 26 . The studies of the isoscalar-isovector ratio 



indicated significant flavor asymmetry when the quarks are light. The difference comes 
from the process with quarks propagating backward in time. In order to understand 
the meaning of the sum-rule violation, we should rely on quark-parton models. 

Proposed theoretical ideas in the 1970's and 1980's are the diquark model [^] and 
the Pauli exclusion effect [|r^, |TI], ^ which were originally intended to explain the old 
SLAG data. In the diquark model, the violation is expected due to the vector-diquark 



admixture. Even though earlier results P, 28 1 seemed to be in agreement with the 



SLAG and NMC data, the deviation from the Gottfried sum becomes very small if the 
virtual-photon interaction with a quark inside the diquark is taken into account with a 
realistic mixing factor between vector and scalar diquarks ^9) . According to the Pauli 



blocking model, uu pair creations are more suppressed than dd creations because of 
the valence u quark excess over valence d in the proton. However, the effect would 
not be large enough to explain the NMC result because a naive counting estimate is 
u/d = 4/5. Furthermore, it was found recently |27] that antisymmetrization between 
quarks could change the situation. If its effects are combined with the Pauli-exclusion 
ones, the u distribution could be larger than the u. 

On the other hand, mesonic models seem to be the most popular idea for explaining 
the NMC result and the flavor asymmetry, at least by judging from number of publica- 
tions. Because of the difference between u and d in virtual pion clouds in the nucleon, 
we have the flavor asymmetry | 



m 



or Tc^p. Because the vr 



|. For example, the proton decays into vr+n 
has a valence d quark, these processes produce an excess of 
d over u in the proton. This model was further developed by including many virtual 
states |3^. Combined mesonic and nuclear-shadowing effects were studied in Ref. [0. 
In the early stage of these models, about a half of the NMC violation was explained 
by the virtual states. In the Adelaide model [^], the NMC deficit was explained by 
adding the Pauli exclusion effect. On the other hand, there is a possibility of explaining 
the whole violation within the mesonic model by considering different ttNN and vrA^A 
form factors 



3^, or by including many virtual states 



Recently, off-shell pion 



m 



but they did not change the pionic contribution 



effects were studied in Refs. 

to u — d significantly. The mesonic mechanism can be described also in chiral models 

m, I3PI, m. ra 1321, ra, 1^, !§], |7[. in the chlral field theory with quarks, gluons, 

m 



and Goldstone bosons [^, ^, ^ 
photon interaction with the pions. 



the flavor asymmetry comes from the virtual 
The obtained results also indicated a significant 

a fraction 



deviation from the Gottfried sum. In the chiral soliton models [p8|, p9|, H 



of the nucleon isospin is carried by the pions, and the deviation from the sum is given 
by the ratio of moments of inertia for the nucleon and pion. A possible relation to the 
a term was also discussed in the chiral models 



41, H3 



|. The virtual mesons could 
modify not only the x distribution ulx) — d{x) but also Q^ evolution of Iq at relatively 
small Q"^ 



m 



u. 



Although it is usually thought to be very small, isospin-symmetry violation, e.g. 
was studied in Ref. 



7^ dn, 



5C, 51 



In order to distinguish the isospin-symmetry 
violation from the flavor asymmetry, we should investigate neutrino reactions, the Drell- 
Yan p-n asymmetry, and charged-hadron production. On the other hand, shadowing 

H to find nuclear 



56, 57 



effects in the deuteron were investigated ||5^, ^ 
corrections in extracting the neutron F2 from the deuteron data. Although there are 
uncertain factors in nuclear potential, the obtained correction to the sum is about 5Ig = 
—0.02. We should mention that it varies depending on the shadowing model. However, 
the correction is a small negative number (except for the pion excess model). If it is 
taken into account, the NMC deficit is magnified! There are also papers on parton- 



transverse- motion corrections [|59| , |60| . It became possible to make flavor decomposition 
in parametrization of antiquark distributions. With the NMC and NA51 data, new 



parametrizations of parton distributions were studied [61, 32 



64] 



The flavor asymmetry in the nucleon could be related to other observables. Nu- 
clear modification of the u — d was investigated in a parton- recombination model [Q . 
Because of the difference between u and d quark numbers in neutron-excess nuclei, 
uu and dd recombination rates are different. This mechanism produces a finite u — d 
distribution in a nucleus even if it vanishes in the nucleon. On the other hand, a rela- 
tion to spin physics was studied pS, EO, E^ . For example, if the Pauli blocking is the 



right mechanism for producing the asymmetry, it also affects the spin content problem. 
Because u[ is larger than u^ in the quark model, a u\. excess over u\ is expected. This 
could be one of the interpretations of the proton spin problem. 

We introduced various theoretical models. In order to distinguish among these mod- 
els, we need theoretical and experimental efforts, in particular by studying consistency 
with other observables. 

The NMC flavor asymmetry can be checked by other experimental reactions. The 
best possibility is the aforementioned Drell-Yan process. We have already explained 
the existing data. Theoretical analyses of the Drell-Yan p-n asymmetry are discussed in 
Refs. [^, ^ ^ ^ |7T], |72|, |73l. The asymmetry should become much clearer by the 
Fermilab-E866 experiment. Charged-hadron- production data in muon scattering by 
the EMC [0, |7^ were analyzed for finding the u — d |j7^. At that time, experimental 



errors were not small enough to judge whether or not the flavor distributions are 
symmetric. However, the recent HERMES measurements show more clearly the NMC 
type flavor asymmetry [|r^. On the other hand, W^ and Z° production can also be 
[T^ , |78| , [79| . Even though the W production is not very sensitive to the 



used 21 



u/d asymmetry in the p + p reaction, the u/d ratio can be measured in the p + p 
[[f9[| . Quarkonium production is usually dominated by the gluon-gluon fusion process; 
however, the u/d could be measured in the large \xf\ region |^ if experimental data 
are accurate enough. Neutrino scattering is another possibility. Combining neutral- 
current and charged-current structure functions, or combining different ones Fi, F2, 
and F3 for a practical purpose, we could obtain the u — d distribution [^, ^ . 

In the following sections, we summarize theoretical and experimental studies on the 
Gottfried sum rule and on the antiquark flavor asymmetry u — dm the nucleon. Future 
experimental possibilities are also discussed. 



2 Possible violation of the Gottfried sum rule 

First, the Gottfried sum rule is derived in a naive parton model. Earlier experimental 
results by the SLAG, EMG, and BGDMS are explained. Then NMG experimental 
results are discussed. We also comment on recent HERMES data. As an independent 
experimental test of the NMG flavor asymmetry, existing Drell-Yan data are shown. 

2.1 Gottfried sum rule 

The Gottfried sum rule is associated with the difference between the proton and neutron 
F2 structure functions measured in unpolarized electron or muon scattering. Because 
there is no fixed neutron target, the deuteron is usually used for obtaining the neutron 
F2 by subtracting out the proton part with nuclear corrections. 

The cross section of unpolarized electron or muon deep inelastic scattering is cal- 



culated by assuming the one-photon exchange process in Fig. |0 
da = ^ = y y (27r)^ 6\k + p-k'-p^) 






where the matrix element is 



M{ep ^ e'X) = u{k\ A')e7^u(fc, A) ^ < X \ e J,(0) \p,a> . (2.2) 

M and m are the proton and lepton masses, k and k' (A and A') are initial and 
final lepton momenta (helicities), and J^ is the electromagnetic current. The proton 
momentum and spin are denoted by p and a, and p^ is the momentum of the hadron 
final state X. The notation Ylpoi indicates that spin average and summation are taken 
for the initial and final states respectively. From these equations, the cross section is 
expressed by a leptonic current part L'^'^ and a hadronic one W^u'- 

d(J = 7TT7 TT7 -^ W^^ — — , (2.3) 

where a is the fine structure constant, s is given by s = (p + /c)^, £" is the scattered 
lepton energy, and Q^ is defined by Q^ = —q^. Throughout this paper, the convention 
^9m = 9ii = 922 = 5'33 = +1 is used so as to have p^ = Po — P^ = M^- The lepton 
tensor can be calculated as 

L^"' = J2 [u{k',X')Yu{k,X)]* [u{k', X')Yu{k,X)] 

\,X' 

= 2 (PA;'^ + k'^k" - k ■ k'g"'') , (2.4) 



in the unpolarized case. The hadronic part is given by 

X a 

ImT. I ^'^ ^^''^ < P' ^1 MO' ^-(0)] IP' ^ > • (2.5) 



AttM- 



Using hght-cone variables q^ = (g°±g'^)/-\/2 with q = (z/, 0, 0, —\/v'^ + Q^) and v = E— 
E', we have g"*" = —Mx/\/2 =finite and q" = ^J2v — i> oo in the Bjorken scahng hmit, 
Q^ — > C)0 with finite x. The exponential factor becomes e*'''^ = e*'' ^ e*'' ^ . Because the 
q^ part is a rapidly oscillating term, the integral vanishes except for the singular region 
of the integrand according to the Riemann-Lebesque theorem. Therefore, the integral 
is dominated by the light-cone region ^"'" ?a 0. In the deep inelastic lepton scattering, 
we can probe light-cone momentum distributions of internal charged constituents in 
the proton. The formal approach for analyzing the hadron tensor is to use operator 
product expansion. It is discussed in section |^ in explaining QCD corrections to the 
sum rule. Here, we do not step into the details and simply discuss general properties. 
Using parity conservation, time-reversal invariance, symmetry under the exchange of 
the Lorentz indices /x and z/, and current conservation, we can express the hadron tensor 
in term of two structure functions W\ and W^'- 

W,. = -W, (,,. - M^) + ly.^ [p, - ^g,) [p. - ^g.) . (2.6) 



From Eqs. ( ^I3|) , (|2.4|) , (|2.6| ), the cross section becomes 



7 2 



dQdE' 4^2 sin^ ^ 



(2.7) 



2 

Scaling structure functions Fi and F2 are defined in terms of Wi and W2'- 

Fi = MWi , F2 = vW2 . (2.8) 

The Fi is associated with the transverse cross section, and the F2 is with the transverse 
and longitudinal ones. In the Bjorken limit, two structure functions are related by the 
Callan-Gross relation 2xFi = F2. In the parton picture, the deep inelastic process 
can be described by virtual photon interactions with individual quarks with incoherent 
impulse approximation. It is supposed to be valid at large Q^ in the sense that virtual- 
photon-interaction time with a quark is fairly small compared with the interaction time 
among quarks. Then, the leading-order (LO) or DIS-scheme structure function F2 is 
given by quark-momentum distributions in the nucleon: 

F2{x,Q^) = Y,elx[q,{x,Q^) +q,{x,Q^)] , (2.9) 



where i denotes the quark flavor. In the next-to-leading order (NLO) except for the 
DIS scheme case, the gluon distribution also contributes to F2 through the splitting 
g — *■ qq. With the assumption of isospin symmetry in the nucleon, parton distributions 
in the neutron could be related to those in the proton. The d-quark distribution in the 
neutron is equal to the u-quark distribution in the proton [un{x, Q"^) = dp{x, Q"^)] and 
in the similar way for other partons [dn = Up, Un = dp, dn = Up, and etc.]. Hereafter, 
the parton distributions are assumed as those in the proton except for section ^3 



where possible isospin-symmetry breaking is discussed. Then, the difference between 
the proton and neutron structure functions is given by 

F|(x, Q2) _ F^ix, Q') = ^x K(x, Q2) _ d^^x, Q2)] + 1^ [^ix, Q^) - d{x, Q^)] . 

(2.10) 

The valence-quark distributions should satisfy 

/ dxu^{x,Q^) = 2 , / dxdy{x,Q'^) = l , (2.11) 

Jo Jo 

due to the proton and neutron charges, / dx{2uv — dv)/S = 1 and / dx{2dy — Uv)/3 = 
where elastic scattering amplitudes are expressed in the parton model by considering 
an infinite momentum frame. Substituting Eq. ( |2.11D into Eq. ( |2.10D and integrating 
over the variable x, we obtain 

/ -^[F^{x,Q')~F-{x,Q')] = - + - dx[u{x,Q')-d{x,Q')] . (2.12) 

Jo -^ "->•-> Jo 

If the sea is flavor symmetric u = d, the second term vanishes and it becomes the 
Gottfried sum rule |]^: 

f^ dr 1 

/ -[Fi{x,Q')-F,-{x,Q')] = - . (2.13) 

Jo X 

As it is obvious in the above derivation in a naive parton model, there is a serious 
assumption of the flavor symmetry in the light antiquark distributions. Therefore, it 
is not a rigorous one like the Bjorken sum rule. Even if violation of the sum rule 
is found in experiments, there is virtually no danger in the fundamental theory of 
strong interactions, quantum chromodynamics. It is nevertheless interesting to test it 
because its violation could suggest an SU(2)-flavor asymmetric sea in the nucleon as it 
was found in the neutrino-induced dilepton production in the case of SU(3). Because of 
small u and d quark masses, large u/d asymmetry cannot be expected in perturbative 
QCD. Therefore, a possible sum-rule violation gives an opportunity for learning more 
details on internal structure of the nucleon. 



2.2 Early experimental results 

Because the small x region could have a significant contribution to the sum rule, it was 
not possible to test it until recently. The minimum x is restricted by the lepton-beam 
energy E as min{x) = Q'^/2ME, where Q^ should not be smaller than a few GeV^ 
in order to be deep inelastic scattering. The first test of the sum rule was studied at 
SLAG in the 1970's. The electron-beam energy is 4.5-20 GeV so that the smallest x is 
about 0.02. Targets are hydrogen, deuterium, and heavier ones. The data are taken in 
the X range from 0.02 to 0.82 for the hydrogen and deuterium targets. The Q^ varies 
depending on the x region, but it is from 0.1 GeV^ to 20 GeV^. In the 1975 analysis |P, 
the data with 0.02< x <0.28 are combined with previous data in the extended range 
X < 0.82. The neutron structure function is extracted by taking into account Fermi 
smearing effects: F^/F|' = (S'F^ — -F|')/F|', where S is the Fermi smearing factor. The 
difference becomes F|' — F2 = 2F|' — 5*^2'^. We define a Gottfried integral by 

/ —[F^{x,Q^)-F^{x,Q^)] = lG{Xmin.Xma.) • (2.14) 

J T ■ ^ 

According to the SLAG data in 1975 H, it is 

/g(0.02, 0.82) = 0.200 ±0.040 (in 1975) . (2.15) 

It should be noted that the integral contains various Q^ data ranging from small Q^, 
where perturbative QGD may not be valid. In any case, it is interesting to find a 
significantly smaller value than the Gottfried sum 1/3. Therefore, there was earlier 
indication of the sum-rule violation in the SLAG data. In fact, the Pauli-blocking and 
diquark models were proposed, just after the SLAG finding, for explaining the possible 
deficit in the sum. However, it was not conclusive enough to state that the sum rule is 
violated experimentally due to a possible large contribution from the smaller-x region. 
Next experimental data came from EMG measurements at GERN by deep inelastic 



muon scattering on the hydrogen and deuterium |0. The muon-beam energy is 280 
GeV, and the measured kinematical range is 0.03< x <0.65 and 7< Q^ <170 GeV^. 
The neutron structure function is extracted from the deuteron data by taking into 
account the smearing effects due to the nucleon Fermi motion. The Hulthen and 
Paris wave functions are used in the 1983 and 1987 analyses to estimate the smearing 
correction. The mean Q^ in the data depends on x, and it ranges from 10 GeV^ at 
x=0.03 to 90 (80 in 1983) GeV^ at x=0.65. Using the Q^ averaged data at each x, 
they obtained 

/g(0.03, 0.65) = 0.18 ±0.01 (stat.)± 0.07 (syst.) , (2.16) 

in 1983. The distribution Ff — F^ is extrapolated into the unmeasured regions by 
using a function F^ — F^ = Ax^'^{l — x)"(l ± (3x), where the constants a and (3 are 
obtained from the data. The 1983 EMG result in the whole x is then given by 
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Jg(0,1) =0.24 ±0.02 (stat.)± 0.13 (syst.) (in 1983) . (2.17) 

In the 1987 report, these values became 

/g(0.02, 0.8) = 0.197 ±0.011 (stat.)± 0.083 (syst) , (2.18) 

and 

/g(0,1) = 0.235 +°:J^° (in 1987) . (2.19) 

It should be noted that different Q^ data are collected to get the integral, whereas the 
sum rule is valid at certain Q^. The above result could be consistent with the sum 1/3 
within the experimental error; however, it is also smaller as the SLAC data indicated. 
Another muon group at CERN, BCDMS, also obtained the sum by analyzing muon 
scattering on the hydrogen and deuterium |T^. The muon-beam energies are 120, 200, 



and 280 GeV. The kinematical range is 0.06< x <0.80 and 8< Q^ <260 GeV^. The 
structure function ratio is obtained with the smearing factor calculated with the Paris 
wave function for the deuteron. Integrating the distribution F^ — F^, the BCDMS 
obtained 

/g(0.06, 0.8) = 0.197 ± 0.006 {stat.) ± 0.036 {syst.) (in 1990) , (2.20) 

at (5^=20 GeV^. The larger-x(>0.8) contribution is negligible, and the smaller-x(<0.06) 
one varies from 0.07 to 0.22 by considering the behavior F2—F2 oc x" with 0.3< a <0.7. 
Because of the large uncertainty from the small x region, they did not quote the in- 
tegral value in the whole range of x. Due to the possible small-x contribution, they 
concluded that it could be consistent with the sum 1/3. 

2.3 NMC finding and recent progress 

Although the earlier data suggested violation of the Gottfried sum, it was not conclusive 
enough because of the large errors and a possible large contribution from the small x 
region. In the NMC experiment, the kinematical range was extended to the small 
X region. The NMC obtained 90 and 280 GeV muon scattering data on hydrogen 



and deuterium targets at CERN |T^. The kinematical range is 0.004< x <0.8 and 
0.4< Q^ <190 GeV^. The difference of the structure functions is calculated by 

FP-F^ = 2F^ '[Z , (2.21) 

where the ratio F2 / F2= 2F2/F2 — 1 is determined by the NMC experiment, and 
the absolute value of the deuteron structure function F^ is given by a fit to various 
experimental data. Nuclear corrections such as the Fermi motion in section |2.2| are not 
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Figure 2.2: Experimental history of the Got- 
tfried sum rule. The SLAG and BGDMS in- 
tegrals are evaluated in the region of 0.02< 
X <0.82 and 0.06< x <0.8 respectively. Be- 
cause unmeasured regions are not included, 
they are shown by the open circles. 
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The NMC reanalyzed the integral by using a new parametrization for F2 including 
their own data and revised F^ jF^ ratios. Their result in 1994 is 

/g(0.004, 0.8) = 0.221 ± 0.008 (stat.) ± 0.019 {syst) (2.25) 

at Q2=4 GeV^. The larger x contribution becomes /g(0.8, 1.0) = 0.001 ± 0.001. The 
smaller x one is /g(0, 0.004) = 0.013 ± 0.005 by the extrapolation Fl - F^ = ax'' with 
a = 0.20 ± 0.03 and b = 0.59 ± 0.06. Then, the overall integral is 

/g(0, 1) = 0.235 ± 0.026 (in 1994) . (2.26) 

The sum is consistent with the previous NMC result; however, the error is slightly 
larger due to more extensive examination of the systematic uncertainties. 

In the HERMES experiment [0, the positron beam energy is 27.5 GeV and hydro- 



gen, deuterium, and ^He gas targets are used. The ratio F|/F^ is extracted from the 
unpolarized hydrogen and deuterium data. The measured kinematical range is 0.015 < 
X < 0.55 (averaged in each bin) and 0.4 < Q^ < 11 GeV^. Because the obtained ratios 
agree with the NMC results, the HERMES experiment seems to support the sum-rule 
violation. However, the sum Iq is not reported yet. On the other hand, a clearer indica- 
tion of the flavor asymmetry is given in semi-inclusive data. As it is discussed in section 



5^ , the charged-hadron production ratio r{x,z) = (N^^ — N^'^ )/{NP'^ — N^'^ ) is 
also related to the u/d asymmetry. The data analysis |T^ clearly favors the NMC 
expectation rather than the flavor symmetric one. 

Because of the small errors, the NMC 1991 result is the first one which made 
us realize that the Gottfried sum rule is actually violated. It strongly suggests the 
flavor asymmetry in the light antiquark distributions, namely a d excess over u in the 
proton. After the NMC finding, many theoretical papers are written on this topic and 
independent Drell-Yan experiments are proposed at CERN and Fermilab. Some Drell- 
Yan experimental results were already taken at Fermilab and CERN. They indicate 
also the NMC type flavor asymmetry. These Drell-Yan results are discussed in section 



2.4 Small x contribution 

One of the reasons why the Gottfried sum rule was not investigated in detail in the 
1970's and 1980's is the lack of small x data, which may contribute significantly. The 
smallest x point of the NMC data is 0.004. Their analysis indicates that the small 
X contribution is Jg(0, 0.004) = 0.013 ± 0.005, which is merely 4% of the sum 1/3. 
In evaluating the integral, they extrapolate the data by using the fitting F|' — F^ = 
0.20x°'^^ to the experimental data. However, it is not very obvious whether the small 
X contribution is so small. Slight variations of the NMC small x data could make a 
significant change in the integral as it is obvious in Fig. BTB. 
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The small x contribution was investi- 
gated in Refs. [^ and ^. Three MRS- 
group (Martin-Roberts-Stirling) parametriza- 
tions, which were available in 1990, were stud- 
ied [ig. They are HMRS-B, KMRS-BO, and 
KMRS-B_ which are fit to various experimen- 
tal data without any small x constraint for the 
HMRS-B, x^ type sea-quark and gluon distri- 
butions in the limit x —>■ for the KMRS-BO, 
and x~°'^ for the KMRS-B_. Comparison of 
these parametrizations with the NMC exper- 
imental data is shown in Fig. 2^ [21]. It 
indicates that the parametrization curves are 

consistent with the data points, and yet they satisfy the Gottfried sum rule without 
any flavor asymmetry for the u and d quarks. The NMC raw data I g {0.004, 0.8) is 
shown by the filled circle with an error bar. It is also consistent with the three predic- 
tions. The only difference comes from small x behavior of F|' — F^. According to the 
KMRS-BO, valence-quark distributions at small x behave like x{uy + d^) ~ x 
xd„ ~ x°'^^. The small x fall-off x°'^^ is much slower than the NMC one x 



Figure 2.3: 
fromRef. ^ 



Small X contributions (taken 



0.27 



0.59 



and 
which 

makes a significant contribution from the small x region. In fact, three parametriza- 
tions have /g(0, 0.004)=0.07— 0.11 so that the missing 10% strength could come from 
the smaller-x region. More recent parametrizations are discussed in section |2.6| . 

Therefore, it is not definite whether the small x contribution is relatively small 
as suggested by the NMC. In the HMRS-E case, the sum 1/3 can be reached if the 
integral region is extended to very small x ~ 10^^°. This is an unrealistic number for 
experimental measurement. However, as it is obvious from Fig. 2 of the KL paper [ ^ , 
the small-x contribution should become obvious at a; ~ 10~^. There is an experimental 
possibility of measuring F2 at such small x by accelerating the deuteron at the Hadron- 
Electron Ring Accelerator (HERA) in Hamburg. However, it is not clear whether such 
experiment could be realized at HERA. Therefore, the best way of testing it, at least 
at this stage, is to use other experimental processes. The NA51 experimental data 
1^ support the NMC conclusion, a d excess over u in the nucleon. More complete 



information will come from the Fermilab-E866 experiment in the near future [^ . 



2.5 Nuclear correction: shadowing in the deuteron 

Because there is no fixed target for the neutron, the deuteron is usually used for 
measuring the neutron structure function F2. In the NMC analyses, the deuteron 
and proton structure-function ratios are measured and they are related to the proton- 



neutron ratio by F'r 



VF| 



2F2D/F2 — 1. Together with world- averaged deuteron 



structure functions, the difference F|' — F2 is calculated by Eq. ( |2.21|) . To be precise, 
the NMC result can be compared with the Gottfried sum only if there is no nuclear 
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Figure 2.4: Virtual photon interaction with 
the deuteron (a) in the impulse approximation 
and (b) in the double scattering case. 

2cr^Af + SayD, where 



The shadowing is traditionally de- 
scribed by the VMD model in particular at 
small Q^. Estimates of its effects on the 
Gottfried sum are found 



57 



m Refs. 

The virtual photon transforms into vector- 
meson states (f), which then interact with 
the deuteron. The hadron-deuteron cross sec- 
tion is given by an individual nucleon term 
in Fig. |2.4| (a) and a double scattering term 
in Fig. p^(b) in the Glauber theory: ayn = 



14 



5a^ 



a. 



D 



vN 



87r2 



d^kTSnik^ 



(2.27) 



The Soik"^) is the deuteron form factor given by the S and D state wave functions: 



SD{k^) = 
as Sa^D 



J dr[u'^{r) + w'^{r)]jo{kr). Then, the virtual-photon cross section is written 
= J^vi^"^/ fv) ^(^vn/i^ + Q'^/M^y. This equation is expressed in the F2 form: 



SFMx) = — 



5(7, 



■D 



Z^ 



/2(1 + QVM2)2 



f2.28) 



The p, uj, and (j) niesons are included as the vector mesons. The most contribution 
comes from the p meson and it is about 80%. With this shadowing correction, the 
deuteron F2 becomes F2D = F2+F2 +5F2D- Because no nuclear correction is assumed 



in the NMC analysis, namely [F2 — -F^ 



2 \NMC 



2F| 



F2D, the Gottfried sum becomes 



'G 



jNMC 

'g 



+ 



-5F2n{x) 

X 



(2.29) 



In Ref. |5^, the VMD model is investigated further by including qq continuum in 
addition to the vector mesons, p, a;, and 0. The model can explain the NMC shad- 
owing data for various nuclei. Applying the same model to the deuteron, they find 
the shadowing correction from (5/^(0.004, 1)=— 0.039 to —0.017 depending on different 



where the 



nuclear potentials. The results qualitatively agree with those in Ref. ||57[ 
Pomeron and meson exchange contributions are added to the VMD one. 

Because the shadowing effects on Iq are more or less same in all realistic models, 
other descriptions are not explained in detail. There are other studies in the Pomeron 
and meson exchange models. We briefly discuss these ideas in the following. For the 
details of formalism, the reader may read the original papers. Historically, the first 
estimate of shadowing contribution to Iq is discussed by the Pomeron exchange model 
3^, Q . A possible way of describing the high-energy scattering in the diffractive region 



is in terms of Pomeron exchange. The virtual photon transforms into a qq pair which 
then interacts with the deuteron. In the diffractive case, the target is remain intact and 
only vacuum quantum number, namely the Pomeron, could be exchanged between the 
qq pair and the nucleons. In the earlier works, the shadowing correction in this model 
was rather large 6Ig ~ —0.08 [^, Q. However, the Pomeron contribution is reduced 
if more realistic deuteron wave functions are used according to Ref. p^. Next, 



meson- 



exchange corrections were investigated in Refs. [54, 57 1. The studied mesons are vr 



and a in Ref. p^ , and p is also included in Ref. [p7 



The formalism is essentially the 
same with the one in subsection |4.3.1| . If the corrections due to the tt, uj, and a mesons 
were taken into account, the NMC result became Iq = 0.29 ± 0.03 |Q. Therefore, 
meson-exchange contributions reduce the discrepancy between the NMC data and the 
Gottfried sum. 
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The VMD contributions are compared with the Pomeron and meson exchange re- 
sults in Fig. \2,.5[ The Pomeron contribution is of the same order of magnitude with 
the VMD effect at (5^=4 GeV^. Because the meson exchange produces extra sea-quark 
distributions, its effects show antishadowing. This fairly large antishadowing cancels 
much of the shadowing produced by the vector-meson dominance and the Pomeron 
exchange. The shadowing due to the Pomeron is rather small at a; > 0.05 com- 
pared with other contributions, and it becomes comparable only in the small x region, 
X < 0.01. Adding these contributions, we show the total deuteron shadowing in Fig. 
2^ . From these results, we obtain the correction to the NMC analysis at (5^=4 GeV^: 
F2 — F2 = (-Ff — F2)nmc + 5F2D- The correction to the sum ranges from (5/0=— 0.026 
to —0.010 depending on the nuclear potential [0. The description in Ref. ^\ is 
consistent with the Fermilab-E665 data [|16|| on the deuteron shadowing. This fact 



suggests the above result of 5Ig is a correct estimate. 
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0.002 
0.000 
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Figure 2.5: Vector-meson-dominance, 

Pomeron, and meson-exchange contributions 
at Q2=4 GeV^ (taken from Ref. |57| ). 



Figure 2.6: Total shadowing in the deuteron. 
The Bochum, Bonn, and Paris wave functions 



are used (taken from Ref. 57 



In the beginning, the estimated shadowing effects on the Gottfried sum were fairly 
large, 5Ig = —0.08. However, the recent numerical values seem to converge into about 
5Ig = —0.02 although there are still uncertain factors due to the nuclear potential. 
In comparison with the NMC value Iq = 0.235 ± 0.026, it is about 10% effect. The 
shadowing studies do not alter the NMC conclusion. However, it has to be taken into 
account carefully because the shadowing magnifies the deviation from the Gottfried 
sum. 

2.6 Parametrization of antiquark distributions 

There are various factors which affect the NMC finding. Even the failure of the Got- 
tfried sum is not undoubtedly confirmed. So present parametrizations of the flavor 
asymmetry u — d is subject to change depending on future experimental results. We 
introduce several parametrizations in the following, but these should be considered 
as preliminary versions. If independent Fermilab Drell-Yan experiments confirm the 
NMC result and the NA51, they should be taken seriously. 
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Early version of the parametrization was proposed in Ref. [jTOl for explaining the 
SLAG data. The Fermilab-E288 collaboration analyzed its Drell-Yan data and obtained 



a parametrization in 1981 [Q. After the NMC measurement, several parametrizations 
have been proposed. The first one is Ref. [^ . They find that the global parametriza- 
tion MRS-B (1988) overestimates the NMC F^ - F^ data at small x even though it 
works for the neutrino structure function xF^ and the Gross-Llewellyn Smith sum rule. 
The differences between the MRS-B and the NMC data are used for finding the flavor 
asymmetric distribution. In the similar way, the difference from the parametrization 
EHLQl is used for finding u — d m. Ref. [^ . 



New global MRS parametrizations were proposed by including the NMC data. The 
total sea-quark distribution at Qg is parametrized as 



xS = 2x{u + d + s + c) = A5x"^(l-x)''s(l + e5X^/^ + 75x) 



(2.30) 



and each distribution is 



u = 


= 0.2 5(1- 


-S)' 


-A/2 


d = 


= 0.2 5(1- 


-5)- 


f A/2 


s = 


= 0.15(1- 


-S) 


•) 


c = 


= 6S/2 , 







with xA = Aax'^'^{1 



X 



\VSi 



1 + Iax) 



(2.31) 

(2.32) 



w 




The parameters are determined by fitting many experimental data. 

The 1993 version is a good example in 
showing flavor symmetric and asymmetric dis- 
tributions in comparison with the NMC F2 data 
|6^ . Therefore, we explain them in the fol- 
lowing. Three possibilities are studied in the 
parametrization: (1) S (same), flavor symmet- 
ric sea u = d, (2) Do (different), asymmetric sea 
u ^ d, (3) D_, asymmetric sea u ^ d with a sin- 
gular gluon distribution. The parametrizations 
are compared with the NMC data, and they ex- 
plain the data fairly well as shown in Fig. |2.7| . 
The flavor symmetric distribution {Sq) deviates 
from the asymmetric ones {Dq, D_) only at small 
x{< 0.01), where the data do not exist. In the 
NMC kinematical region, it is not possible to de- 
tect the differences between these parametriza- 
tions. Because the Sq distribution recovers the 
Gottfried sum 1/3, there is significant contribu- 
tion from the very small x region. However, the sum 1/3 can be reached only at very 



^ NMC (preliminary) 



Figure 2.7: MRS-1993 parametrizations 
are compared with NMC data F2 
(taken from Ref. 



6j) 
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small X ~ 10~^°. The situation should be clarified by the Fermilab Drell-Yan experi- 
ments. The MRS group published new ones after the MRS- 1993 version, in particular 
by including the NA51 asymmetry, HERA data, and the single jet cross sections at 
the Fermilab pp collider. Obtained parameters of the recent 1996 version are listed in 
Table p]l| together with those of the 1993 version. The comparison of the recent one 



with other parametrizations is discussed in the end of this subsection. 



Year 


1993 


1996 


Name 


^0 


Do 


£)_ 


Ri 


R2 


i?3 


i?4 


Q-'o 


4.0 


4.0 


4.0 


1.0 


1.0 


1.0 


1.0 


MS 


215 


215 


215 


241 


344 


241 


344 


As 


1.87 


1.93 


0.054 


0.42 


0.37 


0.92 


0.92 


X 








0.5 


0.14 


0.15 


0.04 


0.04 


Vs 


10.0 


10.0 


6.5 


9.04 


8.27 


9.38 


8.93 


£5 


-2.21 


-2.68 


19.5 


1.11 


1.13 


-1.65 


-2.34 


IS 


6.22 


7.38 


-3.28 


15.5 


14.4 


11.8 


12.0 


Aa 





0.163 


0.144 


0.039 


0.036 


0.040 


0.038 


VA 


/ 


0.45 


0.46 


0.3 


0.3 


0.3 


0.3 


7A 


/ 








64.9 


64.9 


64.9 


64.9 


S 
























Table 2.1: Parameters in the MRS 1993 and 1996 versions. The Ql and A^ are listed in the 
units of GeV'^ and MeV respectively. 



MS 



There are other parametrizations, for example by the CTEQ (Coordinated Theoret- 
ical/Experimental Project on QCD Phenomenology and Tests of the Standard Model) 
group 1^. The recent CTEQ parametrizations included the HERA data, which pro- 
vided information on the small x behavior of the parton distributions. The HERA 
data and others from the CCFR, the Collider Detector at Fermilab (CDF), the NA51 
on u/d are included in the new parametrization analyses. The functional forms of the 
parton distributions are similar to the MRS ones, and they are provided at Q^=2.56 
GeV^. According to the MS version of the CTEQ4 parametrization, the light antiquark 



distributions are obtained as ||63 | 



x{d + u)/2 
x{d — u) 



0.071a; 



^1 
l-x] 



\8.041, 



X^ 
0.501,; ^A8'041q + 30.0a;) 



l + 6.112v^ + x) 



(2.33) 



at Q2=2.56 GeV^. The scale parameter is Ajg =202 MeV. 

In contrast to the above parametrizations, the GRV (Gliick, Reya, and Vogt) model 
supplies input distributions at very small Q^ (~0.3 GeV^). The original motivation 
was to set q{x) = g{x) = at certain small Q^ (= /i^) by allowing only the valence- 
quark distributions. Then, the sea-quark and gluon distributions are considered to 
be produced perturbatively through the evolution from /i^. This attempt is slightly 
modified to the form including valence-like sea-quark and gluon distributions even at /i^ 



D 



MS scheme are plotted in Fig. T^. They are 
the MRS-Rl ||2[, CTEQ4M H, and GRV- 
95 H distributions at Q'^=A GeV^. Because 
the GRV-95 was made so as to agree with 
the older version (MRS- A) of MRS at Q'^=A 
GeV^, the MRS-Rl and GRV-95 distribu- 
tions are almost the same. Even though the 
CTEQ4M agrees with the others in the re- 
gion X > 0.2, it is very different in the small 
X region. The CTEQ4M distribution is much 
smaller than the others. Because the NA51 
Drell-Yan result is taken into account in the 
parametrizations, these distributions are almost the same at x = 0.18. However, we do 
not have enough constraint in the small x region. The only one is the NMC data on 
the Gottfried sum. The small x region should be clarified by future experiments, for 
example by the Fermilab-E866. 

Updated information on the various parametrizations of the parton distributions is 
given at |http://durpdg.dur.ac.uk/HEPDATA/PDl:'] . 



Figure 2.8: x{u- d) distributions at (5^=4 
GeV^ in the MRS-Rl, CTEQ4M, and GRV-95 
parametrizations. 
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3 Expectations in perturbative QCD 

According to the NMC conclusion, the Gottfried sum rule should be violated. In this 
section, we discuss how much corrections are expected in perturbative QCD. First, a 
general treatment of operator product expansion is discussed. Then possible perturba- 
tive QCD corrections to the sum rule are discussed. 

3.1 Operator product expansion 

In order to discuss QCD corrections to the Gottfried sum rule, we introduce operator- 
product expansion which is used in applying perturbative QCD methods to the struc- 
ture functions. The hadron tensor W^i, is expressed as the current product in Eq. 
( |2.5|) . It is known in the light-cone limit ^^ — > that the product is expressed in 
terms of local operators and their coefficients. For example if the current is given by 
JfiiO = 'ipiOlfj.Q'^iO with the charge matrix Q in the free massless Dirac theory, it 
becomes HSlI 



- ( [j^(e/2), j.(-e/2)] + [j.(e/2), -u-^m ) 

"^ j:^^''---^'"s^^^^[9!m^ie)]oi..^^ , (3.1) 

odd n 

where 5(^^) is the 6 function, e{C,^) is a step function: e{C,^) = +1 for ^'^ > and —1 for 
^° < 0, S^aui3 is given by S^au^ =gfiagup + g^lpgva-g^luga|3, and the operator is defined 
by 

o^....M. = ^V^(o)QVi. ■■■i.^(o) • (3.2) 

Virtual forward Compton amplitude is usually analyzed instead of the hadron tensor 
W^y, because it is more convenient to use time-ordered product and to treat interference 
terms. The hadron tensor is related to the imaginary part of the Compton amplitude 
by the optical theorem 2M W^y = (4/7r) ImT^y, where T^^ is given by the time-ordered 
product of currents: 

Here, only the unpolarized case is considered. The amplitude is decomposed into three 
invariant ones [021: 



T^iu{q^, J^) = e^yTiiq^, u) + d^yT2{q^, v) - ie^^ua/3 ^3(?^ ^) , (3.4) 
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where the tensors are defined by e^^ = g^„ - q^iq^/q^ and d^,, = -p^p^q'^/u'^ + iPfj.qu + 
Puqti)/^ — dfiu- The amphtude T^ is the longitudinal one, T2 is the longitudinal plus 
transverse one, and T3 appears only in the weak current case. 

As it is shown in Eq. (|3.1[), a product of current operators could be written by local 
operators and their coefficients. The singular behavior at ^^ ^ can be absorbed into 
the coefficients. Therefore, the Compton amplitude is expanded in terms of possible 
operators. However, infinite number of operators contribute to the amplitude in the 
expansion near the light cone. A convenient way to classify the contributions is to 
introduce twist r, which is defined by the mass dimension of the operator minus its 
spin: T = do — n. For example, the twist for the operator ip'-j^d^^ ■ ■ ■ d^^ip is two 
because the mass dimension of ip is 3/2, the dimension of the derivatives are n, and 
the spin is n+1. In this way, the current product is expanded near the light cone, and 
the amplitude becomes [0, Q 

00 00 

zT(j(OJ(o)) — $^5^Q^(e^/i')r ■ ■ ■ e- ol^...^s^^') , (3.5) 

r=2 n=0 

where C^^^"^,^"^) are called coefficient functions and Oj^^...^^{iJ?) are operators. For 
simplicity, the Lorentz indices /i and z/ are dropped in the above equation. In the case of 
interacting fields, it is necessary to introduce a scale /i^ in renormalizing the operators. 
This is the reason why explicit dependence on the renormalization point /i^ is written 
in the above equation. In this way, the Compton amplitude is factorized into the long 
distance part and the light-cone part which could be handled in perturbative QCD. As 
it is given in Eq. ( |3.1|) , the product of the currents has a singular behavior in the limit 
^'^ — > 0, so that the coefficients could be written in a singular form C^(^^) ~ {1/^'^Y'^/'^. 
Counting dimensions in Eq. ( |3.5|) , we obtain dc = n — dp + 2dj = — r + 2dj where 
do and dj are mass dimensions of the operator and the current. From the dimensional 
counting, we find that the lowest-twist contribution, namely the twist-two, is most 
singular in the operator product expansion. From Eqs. ( p.3|) , ( ^^ , and (|3.5| ), the 
Compton amplitude becomes 

T(g^z.)^5^c:(Q^/.2)o:(/.^)^ , (3.6) 

T,n 

where T(g^, z/) represents T^, T2, or T3. The above C^{Q'^,iJ?) and O.^^iJ?) are defined 
by 

d'i e^^< Cl{i\ f,') e ■ ■ ■ e^" = \Q2/2)n ^« (^'' ^') ' ^3.7) 

- 5^<p,a|0;^...^J/i')|p,(T> = 0;(/i2)p^, ■■■p^„ . (3.8) 
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In relating the Compton amplitudes to the structure functions, the following dispersion 
relation is used: 



T(r,^)-lj_ 



u'du' 



ImT{q^^ v) 



u'du' 

q2/2M V — V 



MW{q'^,u) 



^ X"" In 



dx'x"'-^MWiq^,x') . (3.9) 



Comparing Eq. ( p.6| ) with Eq. ( p.9| ), we obtain moments of the corresponding structure 
function. They are then expressed by the scaling functions: 

dxx''-'F,ix,Q') = J2cl,n{Q^^^')0^^^') , n = 2,4,6--- , (3.10) 

T 

and similar equations for F2 and F3, except that the moments are given by J dxx"'~'^F2 
in the F2 case. Because of the crossing properties of the structure function under 
H *r^ u and x <-> — x, the only even-spin operators contribute in Eq. ( |3.1U| ). The 
moments of the structure functions are thus given by the long-range part, which cannot 
be calculated without resorting to nonperturbative methods such as lattice QCD, and 
the light-cone part which can be evaluated in perturbative QCD. 



There exist only even twists in the expansion Eq. (|3.10|) in the massless quark case. 
Therefore, higher-twist contributions are suppressed by the factor of l/Q"^ compared 
with the twist-two. The Gottfried sum rule is a flavor nonsinglet one. A twist-two 
nonsinglet operator is given by 
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T=2,NS 



,-n-l 



n\ 






D^^ip + permutations 



(3.11; 



where D^ is the covariant derivative D^ = d^ — igT°'A't with eight generators T" of the 
color SU(3) group. 

The renormalization point iJ? is an arbitrary constant, so that physical observable 
should not depend on its scale. This fact leads to a renormalization group equation. It 
can be applied to the coefficients C^^(Q^,/i^) by comparing a renormalization group 
equation for a Green's function with the one for the local operator. In the nonsinglet 
case, it is given by 
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,"2" 



}t,-^K -'"■'" ^^\ 



^NS (Q_ 2 







(3.12) 



where k indicates the structure-function type (fc=l, 2, or 3) and r is omitted for sim- 
plicity. The 'j"''^^ is anomalous dimension of the operator which is related to the renor- 
malization factor of the operator (Z^^ = O^'^^/O^^) by 7"'^^(^) = fiid/d/j,) InZ^^. 
The P function is given by j3{g) = fi{d/dfi)g{fi). The solution of Eq. ( |3.12| ) is 
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(3.13) 
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The anomalous dimension, coefficient function, and j3 function are expanded in a^: 

and (3{g) = —g[Po{g'^/16'K'^) + /3i ((7^/1677^)^ + ■■■]. Then the moments of the structure 
function become 



NS/r^2\ 



Mi:::{Q') = M-^{Qt) 



asm 
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i + c, 
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(3.14) 



where 
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(3.15) 



Because Q^ is an arbitrary scale, it is often convenient to express the above equation 
without Ql: 



MSm = K'n hsiQ' 
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(3.16) 



where A^^ is a constant given by M^^^m = ^^f [1 + Cfc^n^«.(Qo)/(47r)][a.(Qg)]'^". In 
getting various sum rules, ^^^^=1 may be evaluated in the parton model. Then LO and 

NLO anomalous dimensions 7q' and 'j^ are calculated by studying renormalization 
of the nonsinglet operator. In order to obtain B^^^^, we calculate ffist perturbative 
correction to the Compton amplitude and then 0^~ ' (p^ / fJ-'^ , g'^) by considering a 
matrix element of the nonsinglet operator between quark states. From these results. 

Combining these 



the NLO correction to the coefficient function -B^f is obtained 183 



anomalous dimensions and the coefficient, we obtain the NLO correction C^^^^ in Eq. 

(Km. 



3.2 Perturbative correction to the Gottfried sum 

In the previous subsection, it is derived how the moments of a structure function at 
certain Q"^ can be calculated with given moments at Qq by using the prescriptions of the 
operator product expansion and the renormalization-group equation. Before discussing 
the Gottfried sum rule, we ffist check NLO corrections to another nonsinglet quantity, 
for example the Gross-Llewellyn Smith sum rule. It is related to the F3 structure 



functions in neutrino scattering: J dxlF^'^ + F^'^)/2, where F^ = {F^ + F^)/2. In the 



?uN 
3 



parton model without NLO effects, {F^^ + 
so that it's integration over x is three {A^^ 
anomalous dimension vanishes (79 ~ ' =0), 
NLO corrections are given by B^f = —4 and 7" 



F^^)/2 is given by Uy{x,Q 



+ rf^(x,Q2) 

= 3). Because the ffist LO nonsinglet 
the coefficient di becomes di=0. The 

l,NS{-)_ 



fB3l, so that we obtain 



r^NSi- 

•-"3,1 



) 



-4. The notation NS{—) indicates a. q — q type nonsinglet distribution. 
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D 



The NLO correction to the Gross-Llewellyn 
Smith sum rule is about 11%; however, the correc- 
tion to the Gottfried sum is very different. Be- 
cause the NLO term in the coefficient function 
vanishes (-B^f = 0) for the nonsinglet structure 

function ^2^'^, the only contribution is from the Figure 3.1: NLO contribution to the 
NLO anomalous dimension 7;^ . Because the splitting g ^ g. 
structure-function combination in Eq. ( |3.17D is 



given by F!^^ + F^'^ = {u — u) + {d — d) in the leading order, it is a g — g type 
distribution. On the other hand, the Gottfried integrand is given in the parton model 
by (F| — F2)/x = (1/3) [{u + u) — {d + d)] , which is a. q + q type. This difference makes 
the anomalous dimension 7"" ' finite. Even though the LO anomalous dimension 



vanishes in both cases, there is a finite contribution from the NLO process in Fig. pA . 
Namely, the q -^ q splitting becomes possible. Because evolution of the g ± g distribu- 
tions is controlled by the splitting functions Pqq ± Pqg, the g + g evolution is different 
from the g — g one |]8^. Because of baryon number conservation, the first anomalous 
dimension in the NS{—) case has to vanish. However, there is an extra contribution 
from the Pqg (note: Pqq + Pqq = [Pqq — Pqq] + 2Pqq) lu the NS{+) case. The anomalous 
dimension is calculated as 18 



n=l,7VS(+) _ ^ p(l) /„_iN 

= -16{Cl-CFCA/2)PAin = l) 

= -A {Cl -CfCa/2) [13 + 8Ci3)- 271''] . (3.18) 

With the numerical values C(3)=1.2020569..., Cp = (N^ - l)/2N^, Ca = N^, and 
Nc=S, we obtain 7"" ' =+2.5576. In this way, the NLO term becomes Ci = 

^n=i,7V5(+)^^2/?o) = -6(C|-CFCA/2)[13 + 8C(3)-27r2]/(33-2n/). Including the NLO 
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correction, we obtain the Gottfried sum |^, ^M'- 

1 r 

1 + 



2(33 - 2nf) 

0.03552 {nj = 3)\ OsjQ^ 
0.03836 {nj = A)) n 



(13 + 8C(3)-27r2^"^^^'^ 



n 



(3.19) 



The NLO contribution is merely 0.3% at Q^=4 GeV^. It obviously cannot explain the 
large violation found by the NMC The NNLO ct^ correction is estimated recently in 
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(3.20) 



The NNLO correction is about 0.4% at (5^=4 GeV^. We find from these higher-order 
analyses that the perturbative corrections are too small to account the NMC deficit. 

The tiny scaling violation is understood in the following way. The Q^ dependence 
comes from the difference between the flavor-diagonal and nondiagonal splitting pro- 
cesses. Because there are two identical particles in the flavor- diagonal case, they should 
be antisymmetrized. If it could be neglected, the Q^ evolution is flavor symmetric and 
there is no scaling violation in the Gottfried sum. However, the above-mentioned 
antisymmetrization provides the very small scaling violation |[49|| . 

We comment on experimental information about possible Q^ dependence in Ref. 
p5[] . The neutron F2{x,Q'^) is obtained from various proton and deuteron measure- 
ments with nuclear corrections. With the F|' parametrization for explaining the NMC, 
HI, or ZEUS data, the Q^ variation 



Ig{Q') = So [ 1 + ci (as/n) + c^ {ajiif ] 



(3.21) 



is investigated. The obtained parameters averaged over the NMC92, NMC95, and 
HI are Sq = 0.242 ± 0.21, ci = -6.00 ± 0.74, and C2 = 40.4 ± 11.1. The result 
indicates large Q^ dependence which cannot be accounted by the perturbative QCD. 
However, the analysis with the ZUES Ff shows rather different values: 5*0 = 0.383, 
Ci = —12.9, and C2 = 76.2. Therefore, accurate information cannot be obtained at this 
stage. Future HERA measurement of F^ at small x is necessary to find the precise Q^ 
variation. 

The perturbative QCD studies show that perturbative mechanisms cannot account 
for the large violation of the Gottfried sum rule. If the violation is confirmed by further 
experiments, the deficit should come from another source, namely a nonperturbative 
mechanism. 
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4 Theoretical ideas for the sum-rule violation 

The NMC results in 1991 and in 1994 indicate a significant deviation from tlie Gottfried 
sum. We showed that the perturbative mechanisms cannot account for the possible 
violation of the sum rule. It is even not clear whether or not the sum rule is in fact 
violated by considering the small x part. A possible way to answer these problems 
theoretically is to use a nonperturbative approach. Various theoretical ideas have been 
proposed for explaining the deficit in terms of explicit flavor asymmetry u — d ^ 0. 
These ideas are discussed in the following. 



4.1 Lattice QCD 

The most fundamental way to treat nonperturbative physics is to use lattice QCD. 
The following discussions are based on Ref . [^ . The forward Compton amplitude Eq. 



( |3.3|) could be computed by taking the ratio of a four-point function and a two-point 
function: 
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<O^(t-r)O^(0)> 



t-ti,i2>>l/AAf]v 



(4.1) 
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Figure 4.1: Twist-two contributions (taken 
from Rcf. 



where r is given hj t = ti — 12, AMm is the mass difference between the nucleon and 
the first excitation state, and ON{t) is the interpolation field for the nucleon. The 
hadron tensor W^i, is then calculated by the inverse Laplace transformation. 

Euclidean path-integral formalism can 
be used for evaluating the four-point func- 
tion. The leading-twist contributions come 
from the diagrams in Fig. [4.1| . Quark prop- 
agators are involved in the diagrams of Figs. 
^[I](a) and (c), and antiquark ones are in Figs. 
ri](b) and (c). Therefore, antiquark contri- 
butions come from either the connected inser- 
tion in Fig. [4.1| (b) or the disconnected one in 
Fig. ^]l](c). We may call the contribution in Fig. [4.1| (b) from "cloud" antiquarks and 
the one in Fig. |4.1| (c) from "sea" antiquarks, so that an antiquark distribution could 
be written as qi{x) = gf (x) + qf{x). In the same way, a quark distribution is expressed 
as qi{x) = qY {x) + qf{x) + ql{x). If the light-quark masses are equal niu = ma, there 
is no contribution to the flavor asymmetry u — d from the sea graphs in Fig. |4.1| (c). 
Then, the contributions become Iq = 1/3 + (2/3) J^ dx[uc{x) — dc{x)\. 

The hadron tensor has not been calculated directly due to a huge numerical task, so 
that three-point function with one current may be investigated. However, the Gottfried 
sum cannot be calculated because the first moment of Ff — F2 cannot be expressed 
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in terms of the matrix element of a twist-two operator. Therefore, real lattice QCD 
estimate of the Gottfried sum is not available at this stage. Instead, scalar matrix 
elements were studied in Ref. 12611 in order to learn about the cloud contributions to 



the u — d number. The scalar charge / cfx^'i! = J d^k{m/E) X^sK s^k,s + d\ ^dk^s] is 
the sum of quark and antiquark numbers with the weight factor m/E, so that they 
could be a measure of the difference u — d. Because we are interested in the cloud 
antiquarks, only the connected insertion (CI) is discussed in the following. In order to 
reduce the artificial lattice effects, it is better to investigate ratios of matrix elements. 
The ratio of isoscalar and isovector matrix elements for the CI is then approximated 



as 



Rs 



{p\uu\p) — {p\dd\p) 



{p\uu\p) + {p\dd\p) 



l + 2jdx[ 



uAx 



dc{x)] 



CI 3 + 2 f dx[uc{x) + dc{x)] 



(4.2) 



Numerical results are obtained by us- 
ing 16'^x24 lattices with f3=6 and the hop- 
ping parameter fi;=0. 105— 0.154. In Fig. |4.2| , 
the obtained ratios Rs are plotted as a func- 
tion of the quark mass niqa. Because the an- 
tiquark number is positive, the ratio has to 
be smaller than 1/3. The cloud antiquarks 
are suppressed in the heavy-quark case, so 
that the ratio agrees with the valence-quark 



R5O.3 




0.5 



Figure 4.2: Ratio of the isoscalar to isovec- 
tor scalar charge (taken from Ref. [p6|| ). 



expectation 1/3 in Fig. |4.2| . As the quark 
mass decreases, the ratio becomes smaller than 1/3. This decrease should be inter- 
preted as the cloud effects. In order to verify this interpretation, we consider a va- 
lence approximation, which means amputating the backward time hopping. The ratios 
with this approximation are shown by the filled circles and they are 1/3 as expected. 
Next, we discuss comparison with the NMC result. The state with the qq clouds has 
higher energy than one of the valence-quark state, which means that the m/E factor 
is smaller. Therefore, we could estimate the upper bound for the u — d number by 
nu - n^ < [< N\uu - dd\N >doud / < N\uu - dd\N >^aience -l]/2. Their result 
extrapolated to the chiral limit is n^ — nj < —0.12 ± 0.05. It indicates a d number 
excess over the u. The obtained result could be used as a measure of the flavor asym- 
metry although the quantitative comparison with the Gottfried sum is not obvious. 
It is, however, interesting to find that the obtained value is consistent with the NMC 
asymmetry in Eq. ( p.24|) . 

Even though there is no direct estimate of the Gottfried sum in the lattice QCD right 
now, there is an indication of the large light-antiquark flavor asymmetry. It is shown 
that the difference comes from the connected insertion involving quarks propagating 
backward in time by studying the isovector-isoscalar charge ratio. Because the flavor 
asymmetry comes from the cloud antiquarks, physics mechanism behind the above 
results is considered as the Pauli blocking and/or the mesonic effects. 
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4.2 Pauli exclusion principle 

Pauli exclusion model was investigated in Refs. [|Toi [ll| for explaining the SLAG data. 



Because the proton has two valence u quarks and one valence d quark, the uu pair 
creation receives more Pauli exclusion effect than the dd pair creation does. This 
results in the difference between u and d in the nucleon. No qualitative calculation is 
done in Ref. |T^ except for a parametrization based on the above intuition. In order to 
explain the SLAG data (/g=0.27 according to the analysis in Ref. |TD[) the following 
parametrization was proposed 

XM = 0.17(1 -x)^° , xd = 0.17(1 -x)^ . (4.3) 

First, we explain how the flavor asymmetry can be calculated in a model even 
though a realistic four dimensional calculation is not available at this stage. A qual- 



itative calculation on the Pauli blocking effects is discussed in Refs. [^ and ||32|. A 
parton distribution in the nucleon is calculated by |^ 



= \Y.j ^^^Pn-{^-^)p-')\<n\4,+ mp>\' , (4.4) 

V n "■ 

where the subscript c indicates a connected matrix element, '0+ is defined by %l)+ = 
7~7"'"?/^/2, and p'„ is the momentum of the intermediate state. It is the probability of 
removing an antiquark g^ with momentum xp'^, leaving behind a state \n >. The 1+1 
dimensional MIT bag model is used for evaluating the antiquark distribution. However, 
a realistic 3+1 dimensional calculation has not been done yet. We estimate the effect 
by a simple counting estimate. 

In the 1+1 dimension, there are three color states for each flavor. Two of the 
three u-quark ground states and one of the three d-quark states are occupied. It is 
possible to have only one more u-quark in the ground state, but two more d-quarks can 
be accommodated. Therefore, expected sea-quark asymmetry is fairly large: d = 2u 
in the 1+1 dimensional bag picture. In the four dimensional case, there are six states 
(three-color times two-spin states) in the ground state. There are four available ground 
states for u-quarks and five states for d-quarks, so that the asymmetry becomes 

— = - (in a naive counting estimate) . (4.5) 

u 4 



Because there is no valence antiquarks in the bag, Eq. (|4.4| ) indicates that the con- 
tribution comes from a quark being inserted, interacting in the bag, and then being 
removed. Therefore, the d excess is related to the distribution associated with a four- 
quark intermediate state f^^x) 

/ dx[dsea{x) - Usea{x)] = / dxfi{x) = I - P2 , (4.6) 

Jo Jo 
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where P2 is the integral of a distribution associated with a two-quark intermediate 
state. Because the u and d distributions are not calculated in the four dimensional 
model, the Pauli contributions are given by a simple function a;"(l — x)^ in Ref. [0. 
The constant a is chosen to match the small x behavior of used valence distributions, 
and (3=7 is taken so that it contributes only at small x. The overall normalization is 
not determined theoretically at this stage. Roughly speaking, we expect to have I — P2 
in the 25% range because of the naive counting estimate in Eq. ([4.5| ). The obtained x 
dependent results are studied together with pionic effects in the following subsection. 
It is shown that the Pauli blocking effects could produce the excess of d quark over 
u. The naive counting in four dimension indicates d/d=5/4:. Unfortunately, qualitative 
X dependence is not calculated except for the 1+1 dimensional model. It is indicated 
that 10% Pauli effects together with pionic contributions can explain the NMC data 
F2 — F2 fairly well |3^ . However, it was found recently that the conclusion should be 



changed drastically if the antisymmetrization between quarks is considered in addition 



[p7[l . The same -u- valence excess, which suppresses the uu pair creation, also produces 
extra diagrams involved in the uu creation because of the antisymmetrization with the 
extra u. These extra diagrams contribute to a m excess over d. According to Ref. 



2(], if the Pauli-principle and antisymmetrization effects are combined, the u could be 



larger than d, which is in contradiction to the NMC conclusion. 

4.3 Mesonic models 

The Pauli exclusion mechanism produces the flavor asymmetry; however, its effects on 
the sum rule do not seem to be large enough to explain the NMC result according to 
the counting estimate. Furthermore, if they are combined with the antisymmetrization 
contributions, we could have a u excess over d. On the other hand, the meson-cloud 
mechanism is the most successful one in explaining the major part of the NMC flavor 
asymmetry. Because a significant amount of papers are written on this idea, we explain 
the model in detail. We first discuss conventional virtual-meson contributions in section 



4.3. 1[ Second, chiral models in a similar spirit are explained in |4.3.2| . Third, possible 



modification of the Q^ evolution due to the meson emission is discussed in [4.3.3| . 



4.3.1 Meson-cloud contribution 

It is well known that virtual mesons play a very important dynamical role in nucleon 
structure, as they have been studied in the context of cloudy bag or other chiral models. 
The proton decays into vr"*" and neutron, tt" and proton, and other states within the 
time allowed by the uncertainty principle. The virtual pion state is essential for ex- 
plaining many dynamical properties, for example the large A decay width and negative 
square charge radius of the neutron. Therefore, it is important to study whether or 
not the mechanism could produce a flavor asymmetry. It should be noted that pertur- 
bative contributions to the antiquark distributions through the gluon splitting into qq 
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should be much larger than the mesonic ones at reasonably large Q^. However, these 
contributions are supposed to be flavor symmetric, and the asymmetric distributions 
{u + d)/2 — s and u — d could be used for testing the meson mechanism. 

The original idea stems from Ref. [B7 



in 1972, so that the process in Fig. ^^ is 
sometimes called "Sullivan process" . The pro- 
ton splits into a pion and a nucleon, and the 
virtual photon interacts with the pion. Anti- 
quark distributions in the pion contribute to 
the corresponding antiquark distributions in 
the proton. Although the idea is interesting, 
it had not been a very popular topic until re- 
cently. It is partly because experimental data 
are not accurate enough to shed light on the 
mechanism. After the NMC discovery, it is 
shown that the pion-cloud mechanism could 
explain a signiflcant part of the NMC flnding 
pO| , pT| , |32| . This idea is developed further by including many meson and baryon states 
Q and by considering different form factors at the meson-baryon vertices |^ so as 




Figure 4.3: Mesonic contribution to q. 



to explain the whole NMC asymmetry. 

The formalism in Ref 
M = vr is derived by replacing the 7* + p 
in section |2.1| by 



87| is the following. The cross section of Fig. |4.3| with 

e' + X formalism 



X vertex in the e + p 
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where -F^atat (t) is the ttNN form factor, g^^j^ is the vrNN coupling constant, and 0* ■ r 
is the isospin factor. The W2 structure function for the pion is deflned in the same way 
with Eq. (|2.6|) by the replacements p -^ p-j^ and M -^ ?7i^. Then, projecting out the 



F2 part, we obtain 
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(4. 



where tm is the maximum energy transfer: tm = ~'"^Ary^/(l — v)- This equation is 
understood by the convolution of the pion structure function with a light-cone momen- 
tum distribution of the pion. The formalism is used for antiquark distributions in the 
same manner. 

The studies of the pionic mechanism used to be somewhat confusing. The direct 
interaction of the photon with mesons being present in the cloud of a nucleon does not 
contribute to the Gottfried sum. This does not mean that the mesons do not contribute 
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to the Gottfried sum and the u — d distribution as explained in this subsection. In 
deahng with this issue, there are two types of descriptions. One is to calculate only 
mesonic contributions to the u — d distribution pOl, ^, E^ and another is to include 



recoiling baryon interaction with the virtual photon in addition [R^, p3|. Both are 



essentially the same. The details of the compatibility are discussed in the following. 

[I. Models with only meson contributions] 

First, we discuss the former approach with the only meson contributions. The 
Seattle ||3^ and Indiana [^ papers proposed pionic ideas for the sum-rule violation 



in this description. The only major difference is the inclusion of p — > vrA processes in 



Ref. 1 31] in addition to the ttNN ones. Relative magnitude and sign of the uNN and 
ttA^A contributions can be understood in the following way. We consider the processes 
p ^ 71^ + n, TT^ + p, TT^ + A^, tt" + A+, and vr^ + A++, where the virtual photon 
interacts with the pions. Assuming the flavor symmetry in the pion sea, we have the 
u — d distributions in the pions: 

{u - d)^+ = -K , {u - 3)^,0 = , [u- d)^- = +K , (4.9) 

where Kr is the valence-quark distribution in the pions. The flavor symmetry assump- 
tion in the pions does not alter our conclusion unless at very small x with the following 
reason. For example, let us consider the u — d distribution at x=0.1. The light-cone 
momentum distribution of the pion is peaked at y ~0.25; therefore, the most important 
kinematical region for {u — (i)^ is at x/y ~0.4. The valence distribution still dominates 
in this region, so that the sea asymmetry in the pion does not matter. Including isospin 
coefficients at the vrNN and vrNA vertices, 

|0;+-r|2= 2 , ICo •rr= 1 , 

\K+-f\^ = \. ICo-f|2 = ^ , |C_-f|2=l, (4.10) 

we have the isospin times the {u — d).,^ factors as 

^ \4>l-T\^{u-d)^= -2 K for the vrNN process , 

^ |0;-f|^(M- J)^= +-Kr forthevrNA . (4.11) 

o 

In this way, we find that the ttNN contribution to -u — d is negative and is partly 
canceled by a positive contribution from the vrNA. The other important factor is the 
meson-baryon vertex form factor. Because the exact functional form is not known, the 
following monopole, dipole, and exponential forms are usually used: 
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(4.12) 



at the MNB vertex. The different parameters could be related, for example, by A^ = 
0.62Ad = 0.78Ae |3l|]. 

Detailed numerical results are shown in 



ttNA (Aj=O.B0) 



dipole form factor 




. ttNN (Aj=0BO) 



Fig. [4.4| , where pionic contributions from the 
ttNN and vr A^A processes are shown [^ . The 
dipole cutoff parameter A2=0.8 GeV [= A^^ in 
Eq. ( ^121) ] is fixed by fitting the {u + d)/2-s 
experimental data. The dotted curves are 
ttNN and vrA^A contributions. As it is shown 
in the naive discussion, the irNN effect is 
negative and it is canceled by the positive 
ttA^A contribution. The total contribution 
with A2=0.8 GeV is shown by a solid curve 
together with those at A2=1.0 and 1.2 GeV. 
It is noteworthy that the total u — d curve is 
not very sensitive to the cutoff although the 
distribution {u + d)/2 — s does depend much 
on it. Integrating the pionic contribution over 
X, we obtain A/g = 2/3 j dx{u — c?)^ = —0.04, which accounts for about a half the 
discrepancy found by the NMC. It is encouraging that the mesonic model gives a rea- 
sonable value for the magnitude obtained by the NMC. Although we discussed only 
the ttNN and vrNA processes, other processes should be investigated. For example, 
kaon. A, S, and S* are added to vr, N, and A in Ref. |35|. The first method is well 
summarized in Ref. 



Figure 4.4: t:NN and ttA^A contributions 
to u — d (taken from Ref. [31]) 



35[| , so that we quote its results in the following. 



Mesonic contributions to an antiquark distribution in the nucleon are given by the 
convolution of the corresponding antiquark distribution in a meson with the light-cone 
momentum distribution of the meson. The contributions are given by the equation 



MB ^^ 



X 



dyfMBiy)-qMix/y^Q ) 
y 



(4.13) 



vr. 



K) and 
, and a^/B i^ the spin-flavor SU(6) Clebsch-Gordan 



where the summations are taken over combinations of meson states M 

baryon states B = {N, A, A, S, S* 

factors. This equation corresponds to Eq. (|4.8|) in the F2 case. The light-cone momen 
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turn distribution of the virtual meson is 

fM{y) = $^/mb(i/) , (4.14) 



B 



with 



;,„(,). |^,£*^i|~)lf_M,^ . ,4.15) 

The integrand factor X(t, 171^, itlb) is given by 

X(t, m-AT, itib) = —t + {niB — m^Y for i? G 8 

_[{mB + mNy-t] [{niB - m^)'^ - t] 



12 m\jm\ 



for fi G 10 , (4.16) 



depending whether the baryon B is in the baryon octet or in the decuplet. The upper 
hmit of the integral is given by 

tm = m%y--^^ . 4.17 

l-y 

Because the coupling constants are relatively well known, the only factor to be 
paid attention to is the MNB form factors. There are recent studies on whether the 
form factor is hard or soft. Instead of stepping into the details, we summarize briefly 
the historical background and the present situation. A hard form factor with the 
typical monopole cutoff 1.0 < A^ < 1.4 GeV is essential for explaining the deuteron 
D-state admixture and nucleon-nucleon scattering experiments. On the other hand, 
softer ones are obtained in quark models: for example A^ ~ 0.6 GeV in the cloudy- 
bag model [^ 0.7 < A^ < 1.0 GeV in a fiux-tube model ||S^. However, a conflicting 



result came from the studies of the flavor asymmetric distribution {u + d)/2 — s. It 
was originally announced that the cutoff should be much softer, A^ < 0.5 GeV [pO[| , 
which contradicts awfully to the OBEP one. Later analysis with renewed experimental 



data show a slightly larger cutoff A^, ~ 0.6 GeV [31|, which could be consistent with 
those in the quark models. In the Jiilich approach, which is discussed in the following 
paragraphs, the obtained cutoff becomes larger A^ ^ 0.74 GeV (note: monopole cutoff 
is estimated by A^ = 0.62Arf with A^ = 1.2 GeV [^) because more meson and baryons 
are added to vr, A^, A and because the following normalization factor Z is taken into 
account. The probability to find the bare nucleon is reduced by the factor Z due to 
the presence of the meson clouds as explained in detail in the model II. The recent one 



in Ref. |^ without explicit baryon contributions indicates a similar value A^ ~ 0.8 



GeV. Furthermore, it is discussed that kinematical regions, which contribute to the low- 
energy NN scattering and deep inelastic processes, are very different in the form factor. 
The discrepancy between the hard OBEP form factor and the soft one should not be 
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taken seriously. Whatever the outcome may be, it does not change our u — d studies 
significantly if the parameter is fixed so as to explain the {u + d)/2 ~ s distribution. 



[II. Models with meson and baryon contributions] 

Next, we discuss the second approach, which includes recoil-baryon interactions 
with the virtual photon in addition to the meson interactions as shown in Fig. J7S. 
This type of description is studied in the Adelaide paper, [^, the Jiilich [^, and Ref. 
^^. In particular, the Jiilich group developed this model by including many meson 
and baryon states. So far vriV, pN, uN, aN, rjN, vrA, pA, KA, K*A, KT,, K*E, KY*, 
and K*Y* states are included. The pions do not contribute to Fr 
T^NN and nNA contributions are given by 
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(4.18) 



The functions fN^y) and fA^v) are pion light-cone momentum distributions, and Z is 
the valence normalization factor Z = 1/(1 + N.,^ + A^r) with probability of finding a 
pion A^(A)^ = Jq dyfN{A){y)- From these equations, the sum becomes Iq = {Z/3){1 — 
iV^/3 + 5A^/3). According to this equation, the failure of the sum rule is not due to 
the photon interaction with the virtual pion but it is due to the interaction with the 
recoil baryons. This may seem contradictory to the conclusion in the first approach. 



However, it is not a paradox as explained in Refs. [33, 




Figure 4.5: Baryon terms in (b) arc included 
in addition to the meson contributions in (a). 
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Table 4.1: Coefficient Ai and Bi in two dif- 
ferent descriptions 
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Adding other contributions from light meson and baryon states, we write the pre- 
vious equation as 



1 f^ - 

Ig = -{Z + y2Ai) , with Ai= dx {ui + Ui - di - di] 

■^ .• Jo 



Sull 



(4.19) 



where the meson and baryon contributions Ai are given in Table |4.1| . It should be noted 
that all the coefficients in the table should be multiplied by the probabilities of finding 
the meson-baryon states in the nucleon. The nucleon "core" satisfies the valence sum 
J dx{uy — dy)core = 1/3 but its probability is reduced by the normalization factor Z due 
to the virtual MB states. On the other hand, the sum could be written in a different 
form. Whatever the normalization mechanism is, the valence sum should be exactly 
satisfied. Therefore, a part of the meson and baryon contributions can be included 
into the valence sum 1/3. Then, the deviation from 1/3 is identified with the flavor 
asymmetry due to the Sullivan processes / dx{u — d) = f dx{u — d)suii 0) 
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1 f^ - 

- (1 + ^ Bi) , with B,= dx {ui - di)suii ■ (4.20) 



Equation ( |4.2(]| ) corresponds to the first approach without the baryon contributions and 
Eq. ([4.19|) to the second method by the Adelaide and Jiilich. As the coefficients are 
listed in Table ^.1| , there is no contribution from the pion and rho mesons to the sum 
in the second approach. Therefore, the violation comes from the normalization factor 
Z and the baryon contributions. Because the distribution u + u — d — d vanishes for 
example in the pion, this is a natural consequence. However, the virtual ttB contributes 
to the renormalization of the valence-quark distributions. Therefore, we may take out 
the pionic renormalization contributions and put them into the obvious valence-sum 
factor 1/3 in Eq. (|4.2CI| ) [l3ll . Then, it becomes apparent that the pion contributes 



to the deviation from the Gottfried sum as indicated in Eq. ([4.2CI|) . Because of the 
flavor symmetry assumption in the MB, the pion and rho are the only light hadrons 
which contribute to the violation. In this way, the two different mesonic descriptions 
are equivalent, and both numerical results have to be the same. 

In the beginning, it was shown that a significant part, approximately a half, of the 

NMC deficit could be explained by the virtual pions. The Adelaide group tried to 

interpret the whole deficit by adding the Pauli exclusion effects. In the Jiilich model, 

it is explained only by the meson model with the additional meson and baryon states 

3^. The other possibility for explaining the whole result is to consider different form 



factors in the nNN and ttA^A vertices [^. It is also discussed in Ref. ||3^ that the 
normalization factor does not affect the meson part because the bare coupling g^ has 
to be used in the wave function of a physical particle in terms of its constituents. Al- 
though there are slight differences among the above models, the meson-cloud approach 
is successful in explaining at least the major part of the NMC flavor asymmetry. 
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According to the above conclusion, we 
should have the antiquark distributions at 
least in the virtual pion and rho for a real- 
istic evaluation of the violation. So far, we 
have been using the distributions in the pion 
measured in the Drell-Yan processes, namely 
those in the on-shell pion. If the distributions Figure 4.6: Off-shell pion structure function 
in the off-shell pion are significantly modified, in the NJL model. 
numerical results in this subsection should be 

reanalyzed. There is an experimental proposal for measuring virtual pion structure 
functions by detecting a recoil nucleon at HERA ^1|]. Here, we discuss model esti- 
mates of the off-shell effects |36, 37]. Tokyo-Metropolitan-University (TMU) [^ and 



Brooklyn [^ groups studied this topic within the Nambu-Jona-Lasinio (NJL) model. 
While the TMU group calculated distributions at a hadronic scale (Q^=0.25 GeV^) 
and they were evolved to larger Q^, the Brooklyn group calculated them at large Q^ 
directly. The Compton amplitude in Fig. [4.6| is calculated with the Trgg vertex given 
by the NJL model. Then the pion structure function is projected out from the ampli- 
tude, and the sea-quark distributions in the proton are calculated in the pion model. 
Although the ttNN cutoff depends much on the off-shell nature of the pion, the ob- 



tained contribution to Iq is not very significant. According to Ref. ||S^, the deviation 
AIg = —0.0557 with the on-shell pion becomes AIq = —0.0586 with the NJL off-shell 
pion structure function. The difference is merely 5 % effect (0.9% in the Gottfried sum) 
because of the cancellation between the irNN and vrA^A in Eq. (|4.11|) . In any case, if 
model parameters are fixed by fitting other distributions such as {u + d)/2 — s, present 
mesonic contributions are not significantly changed because of the off-shell nature. 

4.3.2 Chiral models 

In the previous subsection, we find that the mesonic contributions could explain the 
major part of the Gottfried-sum-rule violation. The difference between vr"*" and tt~ 
production in the process p -^ Bit gives rise to the antiquark flavor asymmetry. The 
pion production ratio it~^:7t^:7t~ is 2:1:0 in the processes p —>■ Ntc [Eq. ( |4.10| )]. However, 
as it is obvious from Eq. ([4.11|) , the contribution is partly canceled by the p —>■ Avr 
process. In order to have a better estimate, other resonances have to be included. Their 
contributions could be included in a more microscopic approach with effective chiral 
models. In such models, the pion ratio 7r+:7r°:7r~ becomes 4:3:2 if they are produced in 
the process g — >■ gvr. We explain this kind of approaches in this subsection. Although 
chiral quark-meson models were studied slightly earlier ||38|, ^, ^, we discuss first a 
chiral-field-theory approach in Ref. |^, ^, ^, ^ because of similarity in its formalism 



to those in the previous subsection. Later, the chiral-meson models are discussed. 

In describing hadron properties at low energies, it is important to explain sponta- 
neous chiral symmetry breaking. As an effective model for describing such a property. 



36 



we have the chiral field theory. This model is used for evaluating the Gottfried-sum 
deficit [^ . Appropriate degrees of freedom in describing low-energy hadron structure 
are quarks, gluons, and Goldstone bosons. The effective interaction Lagrangian is 

C = ^itD^ + V^)Y^ + igJpA^Y7,tlj + --- , (4.21) 

where ip is the quark field and D^ is the covariant derivative. The vector and axial- 
vector currents are expressed by Goldstone-boson fields 

a'^=1(^^^^^^^^^) ' (4-22) 

^7r7V2 + r//V6 n+ 

^ = exp{m/f) , n=^| TT- -ny^/2 + T]/^/6 

K- If 



Expanding the currents in power of H//, we have V^ = 0(11/ f)'^ and A^ = id^U/ f 
0(11/ f)'^. Then the quark-boson interaction becomes Cuq = —{Qa/ f)'^9^J.^1^15'4'■ 

We give an idea how the deficit arises 
in this model by considering the splitting pro- 5 u.d 

cesses in Fig. ^7f[ A valence u quark splits /v^^ f y^''^ 

into TT"*" and d, and subsequently into d, u, ^ r „ .r 

and d in the left figure. It also splits into 7r° u — > — i— ^^ — d u — ». — i — * — u 

and M, then into m, m, and u or into m, d, and 





d. Noting isospin factors and assigning the Figure 4.7: Valence u quark splitting, 
factor a for the splitting probability u -^ dn^, we have the final state 

u -^ a7r+ + ad + -tt" + -u = — u H d + -u-\ d . (4.24) 

2 2 4 4 4 4 

In the same way, the d quark splitting becomes 

_ a Q a 5a 7a 5a_ a - 

a — i> avr + au -\ — vr H — d = —-u -\ — -d-\ — -u + -d . (4.25) 

2 2 4 4 4 4 

In this simple picture, the sum deficit is estimated by taking the difference between u 
and d in the above equations: 

2, -, 2a , , 

AlG = -iu-d) = -— . (4.26) 

The probability a is calculated in the chiral field theory. With the interaction La- 
grangian Cjjq, the splitting function is given by 
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m'^z 



m?iz/{l — z) and the ultraviolet cutoff is taken a chiral-symmetry- 



where t^.. 

breaking scale: A ?«1169 MeV. This equation is analogous to Eq. ( 4.15 ) in the previous 
meson model. Integrating the function over z, we obtain the probability a for the 
u — ^ n^d process: 



87r2/^ 



dze{X^-T{z))z{\n 



A^ + Ml 
riz) + Ml 



+ Ml 



A2 + M2 



T[Z 



+ MI 



(4.28) 



where t{z) = m'^z'^/^l — z) and 0{x) is a cutoff function defined by 6{x) = 1 for 
X > and for x < 0. With the cutoff A =1169 MeV, the probability becomes 
a = 0.083 which leads to Jg = (1 — 2a)/3 = 0.278. If a larger cutoff is taken, for 
example A =1800 MeV, the sum becomes smaller (10=0.252). The obtained deficit 
is qualitatively in agreement with the meson-cloud models in the previous subsection. 



The X distribution is calculated by the convolution |40, 43 



SjAk Sjk6ii] / — / — q^^^\x/{yz))Puk. 



AN) 



qTAv) 



(4.29) 



where the flavor summation is taken for the indices k and /. However, the calculated 
u — d distribution is concentrated in the small x region at the NA51 scale Q^ = 27 
GeV^, so that the model only with the pion has difficulty in explaining the large NA51 
flavor asymmetry at x = 0.18 |Q. The model is compared with the meson-cloud 
results and it is extended to study the strange quark distribution ^^. On the other 
hand, model consistency is studied among different quantities: the asymmetry u — d, 
the s distribution, and quark polarizations [0. 

As another effective model to describe low-energy properties of hadrons, a chiral 
quark-meson model was proposed. This is an extension of the linear-sigma model with 
replacement of the nucleon fleld by the quark fleld. The Lagrangian density is given 
by 



C{x) = ip{x){i-]' ■ d + glci^x) + IT ■ (f){x)'~^^]}ip{x) 
+ ]^d^a{x)d>'a{x) + ^9^0(x) ■ 9^0(x) 



(4.30) 



-^H^? + 



[x] — V 



2l2 



- F„mla{x) 



The meson fields are treated as classical mean fields, in which the quarks form bound 
states. In this model, the nucleon consists of valence quarks and a coherent superposi- 
tion of mesons, and it is generated from mean-field hedgehog solution. However, it is 
known that the hedgehog states are not eigenstates of spin nor isospin. The nucleon 
with definite spin and isospin should be obtained by a semi-classical cranking method. 
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For slow rotations, cranked meson spin and isospin are linear in the angular velocity. 
The moment of inertia is given by valence-quark and pion contributions: X = Xg +Xjr. 
For discussing the flavor asymmetry in the pion model, an important factor is the num- 
ber difference between vr"'" and tt~ in the proton. The difference is equal to the fraction 
of the proton electric charge carried by the pions: Nt^+ — N,^- = ^ (-^3)„- Then, the 
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(4.31) 



The deviation from 1/3 is related to the fact that a fraction of the nucleon isospin 
is carried by the pions. The fraction is expressed by the moments of inertia for the 
nucleon and pion. 

More rigorous derivation of the sum-rule violation, which is similar to Eq. ([4.31|) , 
is given in Ref. [^ where a similar chiral-quark-soliton model is used. The model 
consists of quark and pion fields with the following functional 



T>'nT>il)T>ilrexp 



i d^xi) {i^ - MW" -m)ip 



(4.32) 



with U^Hx) 



oij5T--rv(x)/fTr 



The nucleon is treated in the same way with the previous 



chiral quark-meson model. However, the matrix element Iq = (1/3) < p | O | p > is 
calculated with a plausible operator O = J d^x['ijj\_{x)T3'ijj+{x) — ^jj_{x)T^ip^{x)], where 
-0+ and ip- cire positive and negative energy parts of ip. As a result, the integral 
becomes a similar equation to Eq. (|4.31|) : 
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with the moments of inertia 



X 



A^. 



val 



E 



< JTal m >< m |{r3, P_}+| > 

Em. — En 



N, 



^vp = -^^f {Em, En] A) < n\T3\m >< m|{r3,P_}+|n > 



(4.34) 



The state | m > is an eigenstate of the single-quark Dirac equation and f{Em, E^, A) 
is a cutoff function. The lyai and X„p are valence quark and vacuum polarization parts 
respectively, and the P_ is the projection operator of the negative-energy eigenstates. 
Numerical results depend on the choice of the dynamical quark mass M; the sum ranges 
from /g=0.235 for M=450 MeV to Jg=0.288 for M=350 MeV. Similar calculation is 
discussed in Ref. |^^ by using the NJL model. From these results, we find that the 
chiral models give similar qualitative results to those of the mesonic models in section 
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4.3. 1[ However, we have to be careful in comparing the effective-model results with 
the NMC value. The effective models are considered to be valid at small Q^, typically 
Q^ ~ ^^CD- O'^ ^^6 other hand, the NMC sum is given at (5^=4 GeV^ where the 
perturbative QCD is usually applied. It is not obvious whether these different Q^ 
results could be compared directly. 

We also comment on other studies in the chiral models. In Ref. |^l|, the Gottfried 
sum is assumed to be related to the matrix element: Iq = (1/3) < p\uu — dd\p >. 
Then it is estimated in a soliton model. The obtained result is a simple relation: 
Ig ^ (1/3) (M„ - Mp)/{md - TTiu) il],j8|. In Ref. ||, the m - d is related to the a 
term by a = 23{'mu + mrf)(l/2 + u — d). Choosing the mass rriu + rrid to fit the pion 
mass and the a term, we find u — d = —0.134. There is also an attempt to relate the 
Gottfried sum to kaon-nucleon scattering cross sections [^]. On the other hand, there 
is an instanton model approach ^^. The instanton induced quark-nucleon interaction 
is described by a Lagrangian with terms which do not vanish only for different quark 
flavors. This feature could be related to the observed u/d asymmetry. Because we do 
not discuss the details of these works, the interested reader should look at the original 
papers. 



4.3.3 Anomalous Q"^ evolution 

We have learned that the scaling violation phenomena in the structure functions are 
successfully described by the perturbative QCD (pQCD). We beheve that the pQCD 
can be used in the Q^ region, Q^ > a. few GeV^. In the recent years, the experimental 
data improved much and they provide us an opportunity to investigate beyond the 
pQCD in the small Q^ region. We explained that the lattice QCD cannot handle the 
structure functions themselves at this stage. Studies of small Q^ physics are difficult 
and they are inevitably model dependent. It is nevertheless important to study such a 
Q"^ region by some kind of model in order to learn about underlying physics. For this 
reason, we have discussed the low-energy effective models in the previous subsection. 
In the following, we discuss possible meson effects on the scaling violation at small Q^. 

We found that the mesonic mechanism 
is a strong candidate in explaining the fail- 
ure of the Gottfried sum. Because of sim- 
ilarity between the quark splitting q -^ gq 
in Fig. ^3(a) and the one into a quark 
and a meson {q -^ Uq) in Fig. |4.8| (b), 
the meson emission is expected to mod- 
ify the Q^ evolution. The process (a) to- 
gether with other splitting processes gives rise 
to the standard Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) evolution equations. Because the appropriate degrees of 
freedom at large Q^ are the quarks and gluons, the evolution has to be described by 




■>-n 



Figure 4.8: Q"^ evolution due to (a) gluon 
and (b) meson 11 emissions. 
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the interactions among these fundamental particles, namely by the DGLAP equations. 
However, meson degrees of freedom may become important at relatively small Q^ . 
If the meson effects are included in the evolution equations, special care should be 
taken for double counting. We discuss possible mesonic effects on the Q^ evolution as 



discussed in Ref. 49 



If the meson interactions are taken into account in addition, the DGLAP evolution 
equations are modified as 

Q 



dt 
d 



(l^ = ^ Pm ® 9i + 5^ P^^% ® 9i + 5^ P^a ® n^ 



dt 

d 

— n, = X Pn„,, ®qj + Y^ Piia<i, ®(l3 + Yl ^n-n, ® lib , (4.35) 

3 3 b 

where IIq is a distribution function of a meson, and t is defined by t = ln((5^//i^). Since 
the current interest is the nonsinglet evolution, obvious gluon terms are omitted for 
simplicity. The notation ® indicates a convolution integral: f ® 9 = f^ dyf{x/y)g{y). 
The function Pgu is the splitting probability of the meson 11 into a quark and an 
antiquark, and Puq (-Pnn) represents the 11 emission probability from a quark (from 
another meson 11). Light pseudoscalar mesons are taken into account in the above 
equations. Isospin and charge-conjugation invariance suggests P„no = Pdno = Pnno = 
Pjno for neutral pseudoscalar mesons, Pun+ = Pdw- = PdTT+ = Pun- = Pqn and PdTT+ = 
Put,- = Pun+ = Pdw- = for chargcd pions. Defining the distribution q^ = {u + u) — 
{d + d), we obtain a nonsinglet evolution equation 

^g+ = (Q,, + Q,,-)®g+ , (4.36) 

where Q is defined by the difference between flavor-diagonal and nondiagonal splitting 
functions Q = P^ - P^^. They are given by P^^^, = Pg,g^, (i=j) and P^,^ = P,.q. (j^j). 
The meson terms cancel out in Eq. ( [4.36|) ; however, there are mesonic contributions 
to the splitting functions. 

A contribution of the 11 emission to the splitting function Pq^q. is given by [PgiqJ'^ = 
da'^'^/dt. Here, a'^*^ is the total cross section integrated over k± for absorption of a 
virtual photon and emission the meson H. Its anomalous dimension is then calculated 
by the Mellin transformation: 7"^^ = 7"^" = da'^*'^ /dt = da'^*'^ /dt, 7"" = 7'''^ = 

•' 'JV'JV N I JV'''JV'JV 

(d/dt) i a'J^'^ /2 + cr^'^/6 + a'/^^' /3 j , where a^^ is the Nth moment of the cross section. 
In discussing the evolution of the Gottfried sum, we calculate the anomalous dimension 
'yf, which is the first moment of the splitting function for the flavor diagonal minus 
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the one for the nondiagonal. The mesonic contribution is 

d 
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In this way, evolution of the Gottfried sum due to the meson emissions becomes 



(4.37) 



lGiQ') = At{t,to)lGiQl), A+(t,to 



exp 



l4\t) + lar'^\t)-^ant) 



^^r''(to)+^^r'(to)-^^r(to) 



. (4.38) 



Exphcit expressions for the above cross sections are presented in the Appendix of Ref. 
[ ^ . In the large mass limit M -^ oo, the cross section falls off like aj ~ 1/M^, so 
that massive meson contributions are smaller than those of the light mesons (vr, rj, and 

V). 

We comment on the consistency with the DGLAP equations at large Q^ and on 
double counting. First, the cross section has aj ~ 1/Q^ behavior at large Q^. It 
means that the anomalous dimension has higher-twist behavior. As Q^ becomes suffi- 
ciently large, the meson contribution becomes eventually smaller than the logarithmic 
pQCD effect. Therefore, the above formalism is consistent with the DGLAP evolution 
at large Q^. Second, the meson contribution to the anomalous dimension should not 
be added to the pQCD one in order to avoid the double counting. We may just use 
the meson picture below a certain boundary Q^ and the pQCD evolution above it. 



Numerical results are shown in Fig. ^ , 
where the Q^ dependence of the Gottfried 
sum is shown. Parameters in the model are 
dynamical quark mass rriii, derivative and 
pseudoscalar coupling ratio (7^, and vertex 
cutoff parameters A and A for two different 
couplings. The mass m^ is determined by 
the cutoff A (=A) with a normalization con- 
dition of the vertex function. In the range 
of A=0.4— 0.8 GeV, it is consistent with the 
constituent quark mass m^ ^ Mp/2 ^ Mp/3. 
Therefore, the cutoff is chosen 0.5 GeV. The 
solid and dashed curves are the results for 
5^^=0.0 and 1.0 respectively, and the dotted and dot-dashed ones are those for (7^=0.0 
and 1.0 including dynamical quark-mass corrections. As it is discussed in section ^, the 
perturbative QCD predicts a very small Q^ variation. Therefore, the scaling violation 
is mostly controlled by the mesonic contributions. Although the results depend much 




Figure 4.9: Q^ dependence of the Gottfried 
sum (taken from Ref. ||49||) 
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on unknown parameters, it is evident from the figure that the sum rule is violated due 
to the scaling violation caused by the quark-meson interactions. This is an interesting 
result, which suggests another possibility of explaining the sum-rule failure. Further- 
more, X distribution of Ff — F^ is calculated at (5^=4 GeV^ by using input distributions 
at small Q^ without sea-quark distributions |^ . The results agree reasonably well with 
the NMC F2 — F2 data. At this stage, these results should be considered as naive ones. 
We need further refinement of the formalism, including discussions on applicability of 
the modified evolution equations at small Q^. In any case, the model predicts very 
strong Q^ dependence of the Gottfried sum, so that it can be checked in principle by 
future experiments. 

4.4 Diquark model 

Flavor asymmetry u — d in a. diquark model was already noticed in 1976 0. Later 
it was elaborated to compare with the NMC result in Ref. |^8[ ^. According to 
the ordinary quark model, baryons consists of three pointlike quarks with spin-parity 
(1/2)"'" and charges 2/3, —1/3, —1/3. It is successful in explaining gross properties of 
hadrons. However, it was rather difficult to understand the small ratio F^/F^ ~ 0.29 
at large x and missing SU(6) baryon multiplets 20 in 1970's, although these problems 



could be explained within the usual quark- model framework |^. These difficulties are 
understood in a diquark model. This model is based on the observation that diquark 
degrees of freedom are most relevant for some observables. 

In the SU(6) model |l|], two quarks form twenty-one symmetric states and fifteen 
antisymmetric ones: 6x6 = 21 + 15. In combination with the remaining quark, the 
antisymmetric part 15 becomes 15 x 6 = 20 + 70. If the antisymmetric part 15 does 
not couple to the quark, it is possible to explain the missing states 20. The SU(3) 
content of the representation 21 is expressed as 21 = {6} x 3 -(- {3} x 1, where the 
brackets indicate irreducible representation of SU(3) and the factors 3 and 1 are spin 
degrees of freedom of the diquark. There are SU(3)-sextet axial- vector diquarks and 
SU(3)-triplet scalar diquarks. We introduce a mixing angle F between the vector and 
scalar diquark states. The usual SU(6) model is recovered in the limit F = tt/A. The 
proton state in this diquark model is given by 

p,s, = ±- \ = ±^=< V2V±i{ud)u^-2V±i{uu)d^ + V2Voiuu)d± 
2 / vl8 t 

-Vo{ud)u± V2smTTS{ud)u±V2cosT\ , (4.39) 

where Vm{qiq2) and S{qiq2) denote the vector and scalar diquark states consist of gi 
and g2 quarks, and the subscript m is the spin state. 

Quark and diquark contributions to the structure function F2 are given by jS^ 
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where Ds{Q^) and -Di,2,3(Q^) are scalar and vector diquark form factors. The Di, D2, 
and D3 are defined by tensor structure of the virtual-photon coupling to a spin-one 
particle: 



V = < {(2fc + g)"(7'^'^Di(Q2) _ [(fc + g)-(^/^" + P(7-"]D2(Q') 

+fc''(fc + qn2k + g)"D3(Q')} ei,.(Ai)4^(A2) 



(4.41) 



where ei^j,(Ai) and e2,^(A2) are the polarization vectors of initial and final diquarks 
with helicities Ai and A2. In the limit of pointlike diquarks, the form factors are given 
by Ds{0) = 1, -Di(O) = 1, -D2(0) = I + k, and .03(0) = 0, where k is the anomalous 
magnetic moment. Therefore, it is natural to choose Ds{Q'^) = Di{Q'^) = D2{Q'^) = 
Qo/iQo + Q"^) — D{Q^), as expected from a dimensional counting rule, and Z^s = for 
simplicity. From Eqs. (|4.39|) and (^4.40|) , valence quark and diquark contributions to 
the proton F2 are 



F^,=x 



'^\Ux) + Y^Ux)\s\n'V + ^Ux)oos''V + ^f,{x)cos''VD\Q') 



— -fy ix) + --fy(x)\sin^v(l + ^^ 



D\Q' 



. (4.42) 



New distribution functions fg,s,v{x) are introduced in the above equation, and they 
are normalized as J^ dxfq^s,v{x) = 1. In the same way, the neutron F2 is given by 



F^ 



X 



^lUx) + ^^Ux)\sin'r + ^Ux)cos'r + ^f^{x)cos'rD\Q') 



+ {--fy (x) + --fy(x)\ sin" t(i + ^^ 



D\Q' 



. (4.43) 



In addition, the virtual-photon scattering off a quark inside the diquark is considered 
with a diquark-breakup probability 1 — F'^{Q'^). From these equations, the Gottfried 
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sum becomes p9| 



+ ^|l-F2(Q2)|sm2r . (4,44) 

It should be noted that sea-quark contributions are not taken into account. In other 
words, the sea-quark distributions are assumed to be flavor symmetric. If the nucleon 
consists of a scalar diquark and a quark (F = 0), the integral becomes the Gottfried 
sum 1/3. 

In the earlier investigations [^, |2^ without the breakup term, the diquark 



model seemed to account for the deficit of the Gottfried sum. For example, 
7(^=1/3-0.384 sin^F at Q'^=A GeV^ was obtained in Ref. ||2^. However, it means 
that the two quarks in a diquark act as a single extended object, which never breaks 
apart. This is certainly not realistic. The breakup mechanism plays an important role 
in the Gottfried sum. With the distribution fvuS^) = 12x^(1 — x) used in the study 
of polarized structure function gi and with the assumption F{Q'^) = D{Q'^), the sum 
at (5^=4 GeV^ becomes pi 



1 4 
3 + 9 



,c = - + -(0.12) sin^F . (4.45) 



On the contrary to the previous results, the model produces a positive modification 
to the sum. The parameter F could be taken from other observable such as the ratio 
of the axial vector to the vector neutron /3-decay coupling constant g^. Comparing 
g^ = 1 + (2/3)sm(2F) in the diquark model with experimental value (7^=1.261±0.004, 
we obtain siri^T Ri0.04. Then the sum becomes Iq = 1/3 + 0.002, which is a very small 
positive correction to the sum 1/3. 

There are two factors which changed the early result in Ref. |2^. The first one is 
the addition of the breakup term, and it could result in the positive contribution in 
Eq. ( |4.45| ). The second one is the small mixing angle siri^T ?^0.04, which is consistent 



with (yf^. It should be noted in Ref. ||2^ that the NMC result could be explained if 
the mixing were sin^T ^ 0.27! Sensitivity of the obtained Ig on the assumed functions 
is also discussed in Refs. [p8| , 29|. If the function F{Q'^) is different from D{Q'^), the 



deviation may become negative instead of the positive one in Eq. ( [4.45| ). The results 



depend, of course, on the assumed functions for fvnui^)^ F^Q"^), and D^Q"^); however, 
the diquark effects are all suppressed by the small factor sin'^T. In this way, the diquark 
model predicts a very small deviation from the Gottfried sum, so that it cannot explain 
the NMC results. 
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4.5 Isospin symmetry violation 

Isospin symmetry is usually taken for granted in discussing parton distributions in the 
proton and neutron. In fact, it is assumed [m„ = dp, dn = Up, Un = dp, dn = Up, and 
etc.] in deriving Eq. (|2.10|) . Electromagnetic interactions are weak compared with 
strong interactions, so that typical isospin-violation effects are expected to be of the 
order of the fine structure constant a = 1/137. This is in general true, for example, the 
mass difference of the nucleons is (m„ — mp)/'mp=0.14%. Therefore, we cannot believe 
that the NMC result is explained only by the isospin-symmetry violation in antiquark 
distributions. However, it is worth investigating its contributions to the Gottfried sum 
and to various high-energy processes because isospin-violation effects on the parton 
distributions are not known. This topic is discussed in Ref. fSOl . 

What would happen to the Gottfried sum if the isospin symmetry cannot be as- 
sumed? Without using the isospin symmetry, the sum is expressed as 

^G = l + lj dx{ [A{u{x) + c(x)} + {d{x) + s(x)}]^ 

-[A{u{x)+c{x)} + {d{x) + s{x)}]^^) . (4.46) 



If the antiquark distributions are flavor symmetric and if the s and c terms vanish: 
J dx[sp — s„) = and / dx{cp — c„) = 0, it becomes 

Ig = ^ + ^ dx [qp{x) - qn{x)] , (4.47) 

where qp{x) is the light antiquark distribution in the proton [qp = Up = dp] and g„(a;) 
is the one in the neutron. If the isospin-symmetry breaking were the only origin of the 
NMC finding, Eqs. ( [4.46 ) and ( [4.47 ) could suggest that there are more antiquarks in 



the neutron than those in the proton. If the NMC 1991 data in Eq. ( [2.23| ) is identified 
with Eq. ( [4.47] ), we get J dx[qp{x) — qn{x)] = —0.84 ± 0.014. Because the Adler sum 
rule I A = J[F'2^{x) —F2^{x)]dx/2x = 1 and the Gross-Llewellyn Smith sum rule in Eq. 
([3.17[) are independent of the flavor asymmetry and the isospin symmetry violation. 



these mechanisms cannot be distinguished. We have to find other observables. Future 
experimental possibilities are discussed in section [5.6|. 



4.6 Flavor asymmetry u — d in nuclei 

The NMC finding of the flavor asymmetry can be tested by the Drell-Yan experiments. 
There exist Drell-Yan data for various nuclear targets, so that some people use, for 
example, the tungsten data in investigating the flavor asymmetry [|I^. However, we 
have to be careful in comparing the NMC result with the tungsten data because of 



possible nuclear medium effects. In the analysis of Ref. [^, no nuclear correction is 



made except for the overall shadowing correction. If the nuclear modification in the 
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u — d distribution is very large, the Drell-Yan analysis cannot be compared directly 
with the NMC result. Therefore, it is worth while estimating the nuclear effect in order 
to find whether or not such a modification should be taken into account. 

In discussing antiquark distributions in a nucleus, it is essential to describe the 
shadowing phenomena. In the small x region, it is experimentally observed that nuclear 
structure functions per nucleon are smaller than the deuteron's. There are various 



models in describing the shadowing as they are discussed in section 2.5. The interesting 



point is to find whether there is any nuclear mechanism to create extra flavor asymmetry 



whatever the shadowing model is. It could be possible according to Ref . [^ , in which 
the nuclear modification is calculated in a parton-recombination model. The following 
discussion is based on this investigation. 

First, we discuss the nuclear u — d distribution without the nuclear modification. 
If the isospin symmetry could be applied to the parton distributions in the proton and 
neutron, the distribution per nucleon is given by [u{x) — d{x)]A = —e[u{x) — d{x)]proton- 
It is simply the summation of proton and neutron contributions. The neutron-excess 
parameter e is defined by e = {N — Z)/{N + Z). The above equation indicates that 
the flavor asymmetry has to vanish if the antiquark distributions are flavor symmetric 
in the nucleon. However, it is not the case in the recombination model. 

In the parton-recombination picture, partons in different nucleons could interact 
in a nucleus. These interactions become important especially at small x. Parton 
recombination effects on the antiquark distribution qi{x) are given by 

/^qi,A{x) = Wpp Aqi^pp{x) + Wpn Agj,p„(a;) + Wnp Aqi^np{x) + w„„ Agi,„„(a;) , (4.48) 

where w„j„2 is the combination probabilities of the two nucleons rii and ^2. For exam- 
ple, Wpp is the probability of the proton-proton combination {wpp = Z{Z — 1)/[A{A — 
1)]). The distribution Aqi^nin2{x) is the modification of the antiquark distribution with 
flavor i due to a parton interaction in the nucleon nl with a parton in the nucleon n2. 
If the isospin symmetry can be used, the flavor asymmetry becomes 

x[Au{x) — Ad{x)]A = —{wnn — Wpp) X [Au{x) — Ad{x)]pp (4.49) 

= —e X [Au{x) — Ad{x)]pp , 

where [Au{x) — Ad{x)]pp is the asymmetry produced in the proton-proton combination. 

Next, we discuss how the flavor asym- 
metric distribution is created in this model 
for the simplest situation, u — d = in the nu- 
cleon. In this case, many recombinations can- 
cel each other, and the only remaining term 

. ,, „ „ . figure 4.10: Mechanism of creatine the fla- 

is the foUowmg: , . , 

^ vor asymmetry m a nucleus. 

4Kn C^ x"^ -\- X^ 
x[Au{x) - Ad{x)]A = e——x j dx2xu*{x)x2[uv{x2) - dv{x2)]-, fr , (4.50) 




[X + X2) 
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n 

The recombination contributions are 
evaluated for the tungsten yf^VTuo nucleus 
by using the input parton distributions MRS- 
Dq at Q^=4 GeV^. The obtained results are 
shown in Fig. [4.11| , where the solid (dashed) 
curve shows the x[/S.u — AJ]^ distribution of 
the tungsten nucleus with the flavor symmet- 
ric (asymmetric) sea in the nucleon. In the 
{u — (TjN = case, the positive contribution 
at small x can be understood by the pro- 



Figure 4.11: Created flavor asymmetry 
in the tungsten nucleus by a recombination 
model 



cesses in Fig. 
case, the q{x)G 
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In the {u — d)N 7^ 
q process is the dominant 
one kinematically at small x. Its contribu- 
tion to u{x) — d{x) becomes negative due to 
the neutron excess. On the other hand, the 
qG -^ q{x) process becomes kinematically favorable in the medium x region. Because 
it produces q with momentum fraction x, its contribution becomes opposite to the one 
at small x. 

The above results are obtained at Q^=4: GeV^. Because the factor Kq is proportional 
to as{Q^)/Q^, the nuclear flavor asymmetry may seem to be very large at small Q^. 
However, the quark distribution Uy{x) — dv{x) becomes very small in the small x region, 
so that the overall Q^ dependence is not so significant according to Eq. ( [4.50| ). There 
are merely factor-of-two differences between the asymmetric distribution at (5^=4 GeV^ 
and the one at Q^ ~ 1 GeV^. Considering this factor of two, we find that the nuclear 
modification is of the order of 2%-10% compared with the asymmetry u — d suggested 
by the MRS-DO distribution. Therefore, special care should be taken in comparing 
flavor asymmetry data of the nucleon with the nuclear ones. On the other hand, 
because the Drell-Yan experiments on various targets are in progress at Fermilab, the 
nuclear modification oi u — d could be tested experimentally. The studies could provide 
important clues in describing nuclear dynamics in the high-energy region. 
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4.7 Relation to nucleon spin 

We explained various mechanisms of creating the flavor asymmetry. It is natural that 
there is a certain relationship between the light-antiquark flavor asymmetry and the 



nucleon spin issue. We discuss possible relations in the exclusion principle |66[ and in 
the chiral quark model |^0|, ^Tj . 

One of the ideas in explaining the flavor asymmetry is the Pauli blocking model in 
section E^. This idea could be extended to the proton-spin problem according to Ref. 



[ pq] . In a naive quark model, polarized valence-quark distributions are related to the 
matrix elements of axial charges: 

ul = l + F , ui = l-F , 4= ^^^"^ , 4= ^"^^^ , (4.51) 

where F and D are axial parameters, and the current values are F + iI'=1.2573±0.0028 
F/Z)=0.575±0.016 experimentally. From these equations, fractions of the proton spin 
carried by the valence quarks are 

Au, = ul-ui = 2F , Ad, = dl-di = F-D . (4.52) 

The Pauli blocking mechanism in the flavor case was the following. Because there is 
an extra u- valence quark over d-valence, uu pair creations suffer more exclusion effects 
than dd creations. Substituting numerical values of F and D into Eq. ( [4.51| ), we 
obtain mJ=1.46, m|,=0.54, 4=0-33, and 4=0-67. The proton spin is dominated by the 
ul distribution. Because mJ is signiflcantly larger than u^,, the Pauli blocking could be 
applied to the spin case in the similar way. As a rough estimate, the fraction of the 
spin asymmetry created in the exclusion principle is assumed to be the same with the 
one for the flavor asymmetry 

4^ = ^^ ' (4-53) 

Uv — Uv Uy — dv 

where g^ denotes a sea-quark distribution. The NMC result in 1991 indicates ds — 
Mo=0.14 for the flrst moments m= and d^. Therefore, the difference becomes 



u 



T_4 = -0.14 Am„ = -0.28 F . (4.54) 



We also assume that the exclusion mechanism is applied in the same way to the d 
quark: 

4-4 = -0.14 Arf^ = -0.14 (F-L)) . (4.55) 

From these equations, the flrst moment of gi{x) becomes 

/■' 1 

/ g{{x)dx = —{9F-D){l-0.28) = 0.U, (4.56) 

Jo 18 
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which is in fair agreement with polarized experimental data. The Pauli blocking inter- 
pretation of the proton-spin issue is summarized in the following way. Because of the ul 
excess over ul and the dl excess over dl, the u-quark (d-quark) sea is negatively (pos- 
itively) polarized. However, magnitude of the exclusion effect is expected to be larger 
in the u-quark sea because of ul/ul > d\^/d\. In fact, we have A'u=— 0.28F = —0.13 
and A(i = — 0.14(F — D)=+0.05. The large negative polarization in the u-quark sea 
could account for the spin deficit. 

The relation between the quark flavor and spin can be discussed in the meson 



models, for example in the chiral quark model P0| , ^ . As explained in section |4.3.2 , 
the model produces the d excess over u. In Ref. [0, kaons, eta, and eta prime mesons 
are included in addition to the pions. Suppression factors are introduced for heavier 
meson emissions in comparison with the pion case: e for kaons, 5 for eta, and C, for eta 
prime mesons. Then, after the emission of a meson from the initial proton {2u + d), 
the u — d number becomes 



u — d 



K + 5 ^ 



(4.57) 



The proton spin-up state is |p+) = (1/a/6) (2 \u+u+dJ) — |M+M_(i+) — \u^u+d^)). This 
equation suggests the naive quark model prediction Am = m+ — n_ = 4/3, Ac? = 
d+ — d^ = —1/3, and As = with m+ = 5/3, U- = 1/3, d+ = 1/3, and d- = 2/3. 
Next, we discuss corrections to the quark polarization due to the meson emissions. 
They are expressed by the probability P{q+ -^ g_) for the splitting process q+ ^ q_. 
The probability a is assigned for the process m+ -^ n^d^. P {d^ ^ s_) = e^a is 
the probability of a spin-up d quark flipping into a spin-down s quark through the 
emission of K~^. In this way, the probabilities for the u+ splitting processes are P{uj^ -^ 
{ud)QdJ) = a, P{u+ -^ (■us)os_) = e^a, P{u^ -^ (mm)oM-) = [{6 + 2( + 3)/6]'^a, P{u^ -^ 
{dd)oU-) = [{6 + 2C - 3)/6]2a, and P{u+ -^ (ss)on_) = [{6 - 0/3]^a. Substituting 
these equation and similar ones for other states, we obtain 

M_) - ^P{d+ ^M-) 



rf_)-ip(d+^rf_) 



Au =^ [1 - SP] + ip (m_ ^ «+) + ^P {d_ -. 


5 
u+)--Piu+~ 


4 21 + 452 + 8C^ + 12e2 




3 9 




Arf = - ^ [1 - SP] + ip(«_ ^ rf+) + ^P(d_ 


^ d^) - ^P (u, 


1 6 - 52 - 2(2 - 3e2 
~ 3 9 




As =^P (n„ ^ s+) + ^P (rf_ ^ s+) - ^P (n+ 


^s_)-ip(rf 



(4.58) 

The parameter values, e, 6, and (, have to be determined in order to evaluate above 
quantities numerically. We may simply take the same suppression factors for the K 

50 



and T] production processes. Because rj' is heavier, a smaller value may be taken for 
(. For a rough estimate, we assume e = 6 = —2(. These parameter values are fixed 
by Eq. (|4.57| ) so as to explain the NMC u — d, and they are e = 6 = —2( = 0.6 and 



a = 0.15. These are substituted into Eq. ( [4.58| ) to obtain Am=0.87, Ad = —0.41, and 
As = —0.05. The total quark spin content is Au + Ad + As=0.4 which is close to the 
recent measurements. Comparing these with the naive quark model predictions Au = 
4/3, Ad = —1/3, and As = 0, we find that the u quark polarization is significantly 
reduced due to the meson emission mechanism. Although the above calculation is a 
very rough one, the results indicate that the chiral quark model could also explain the 
proton spin problem. In this way, we find that the fiavor asymmetry problem is closely 
connected with the proton spin issue. 

4.8 Comment on effects of quark mass and transverse motion 

Although the parton model is considered to be valid in the Bjorken scaling limit, 
there could be some corrections from finite quark masses and transverse motion. Such 



corrections are estimated in Ref. |£9|, and the obtained result indicates 6Ig = —0.01 
to —0.02 at Q^=4: GeV^. Later, a more careful analysis indicates a slightly larger 
correction SIq = —0.029 to —0.051 at Q^=3 GeV^ @. It is interesting to find that 
the discrepancy between the NMC result and the sum becomes smaller, but it is not 
large enough to explain the NMC deficit. Because the correction becomes smaller: 
Sla = -0.009 to -0.017 even at slightly larger Q^ (^=io GeV^). This kind of simple 
interpretation could be tested by future experiments. For the details of this topics, the 



interested reader may look at the original papers |59, |60| 
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tion process in hadron-hadron collisions A + 
B — > i^i'X, where i is for example the 



muon, as shown in Fig. |5.1| . The Drell-Yan 
experiments have been used in determining 
quark distributions in a hadron, in particu- 
lar sea-quark distributions. Therefore, it is 

T^. r 1 ^ ,, ^r ideal for examining the flavor dependence in 

rigure 5.1: Drell-Yan process. , ,. , . ,,..,. , , 

the hght antiquark distributions, even though 

there is an undetermined i^-factor in the experimental analysis. The /('-factor is the 
ratio between the measured cross section and the leading-order prediction. The theo- 
retical studies of a. and a? corrections revealed that the /^-factor could be rather well 



explained by the higher-order QCD corrections p5[| . 
Its cross section is given by [^ 

da = ^ y V (2nY5(PA + PB-h-k2- Px) 

X \MiAB^trX)\^^^^^ , (5.1) 

where k\ and k2 are i,~ and t^ momenta, spin summation is taken for the final state 
particles, and spin average is taken for the initial hadrons. Because the matrix element 
is given by 

M{AB^i+i-X) = u{k^Ai)eiAk2A2) ,, , ,, <X\eJ,mAB> , (5.2) 

(fci + k2Y 

the cross section is written in terms of lepton and hadron tensors: 

v/[s-(M2+M2)]2-4M2M| (fci + k2y (27r)32fcio (27r)32fc2o ' ^ ' ' 
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described by the quark- ant iquark annihila- 
tion process q + q — > i~^ + i~ . For exam- 
ple, Fig. 1^ indicates that a quark with 
the momentum fraction Xi in the hadron 
A annihilates with an antiquark with X2 in 
the hadron B. Considering the color factor 
3 ■ (1/3)2 = 1/3, we obtain the LO Drell-Yan 
cross section 



Figure 5.2: Dominant contribution to the 
Drell-Yan cross section. 



'd^ = ^ E ^' [ ^^^(^1' ^') ^t (^2, Q') + qfixu Q') qf{x,, Q') ] , (5.6) 

where Q^ is the dimuon mass squared: Q^ = m?^^, and r is given by r = m'^^^/s = Xia;2. 
The rapidity y is defined by dimuon longitudinal momentum P£ and dimuon energy 
E* in the cm. system: y = {l/2)ln[{E* + Pl)/{E* — PI)]. The momentum fractions 
Xi and X2 can be written by these kinematical variables: Xi = y/re^ and X2 = a/tc"^. 
According to Eq. (p. 61), the process can be used for measuring the antiquark dis- 
tributions if the quark distributions in another hadron are known. For finding the 
flavor asymmetry u — d, the difference between p-p and p-n (practically p-d) Drell-Yan 
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cross sections is useful. Considering the rapidity point y=0 and retaining only the 
valence-sea annihilation terms, we have 



a 



pp 



a 



pn 



QTCa 



-u^{x)u{x) + -dy{x)d{x) 

y y 



4 - 1 

- {uy{x)d{x) + dy{x)u{x)} + -{dy{x)u{x) +Uy{x)d{x)] 

y^T L y y 



at y = , 

(5.7) 



for the proton-proton and proton-neutron cross sections. All the above distributions 
are at a; = ^/T because oi y = 0. From these equations, the p-n asymmetry becomes 



A 



(jPP _ ^pn 



DY 



(jpp + aP"^ 

[Auv{x) — dv{x)][u{x) — d{x)\ + [uy{x) — d„(x)][4u(a;) — d{x)] 

[Auy{x) + dy{x)][u{x) + d{x)] + [uy{x) + dy{x)][Au{x) + d{x)] 



at y = 



(5.^ 



This quantity is very sensitive to the u — d distribution. However, antiquark-quark 
annihilation processes also contribute to the above equation at the rapidity point y=0. 
In this sense, it is better to take large xp (= xi — X2) data so that the antiquarks in 
the projectile do not affect the asymmetry 



_ [4n(a:i) - d{xi)\[u{x2) - d{x2)\ 
[Au{xi) + d{xi)][u{x2) + d{x2)\ 



at large Xp 



(5.9) 



The above discussions are, of course, based on the LO cross section. On the other hand, 
the higher-order corrections are rather large as represented by the if- factor. Therefore, 
it is important to investigate whether the if-factor cancels out in the asymmetry Auy 

The Drell-Yan process has been already used for studying the flavor asymmetry. 
There are existing data by the Fermilab-E288, the Fermilab-E772, and the CERN- 
NA51. Furthermore, detailed studies of the Drell-Yan asymmetry are in progress at 
Fermilab by the E866 collaboration. 

The Fermilab-E288 collaboration measured dileptons produced in proton-nucleus 
collisions. Proton-beam energies are 200, 300 and 400 GeV, and targets are beryllium, 
copper and platinum. The dimuon data are taken in the mass region m^^=4— 17 GeV, 
and they are analyzed by using Eq. ( [5.6| ). The isospin symmetry is assumed for parton 
distributions in the proton and neutron. No nuclear correction is made except for the 
Fermi motion correction. We note that the E288 paper was published before the finding 
of the EMC effect [^. Q^ dependent F2 data from electron and muon scattering are 
used together with a fit F^ jF^ = 1.0 — 0.8x and parametrized antiquark distributions. 
The antiquark part is assumed to be Q^ independent, and they are determined by 
fitting their data: d = 0.548(1 - xY'^^ and u = 0.548(1 - a;)"\ Calculated DY cross 
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sections with these u and d distributions are shown by the dashed curve in Fig. p73 , 
and flavor-symmetric ones are shown by the sohd curve. The E288 data favor a d 
d{l - a:)3-48. 
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Figure 5.3: Slope of rapidity distribution at 

y~^ in the E288 experiment (taken from Ref. 

12 ). Flavor symmetric and asymmetric re- 




sults are shown by the solid and dashed curves. 



Figure 5.4: Drell-Yan cross section ratio 
cw/f/S in the E772 experiment (taken from 
Ref. 18|). Open circles are the data without 
shadowing correction. Three theoretical asym- 
metry results are shown by the solid, dashed, 
and dot-dashed curves. 



Next, the Fermilab-E772 collaboration showed their Drell-Yan results on the flavor 



asymmetry in 1992 |T^. The Drell-Yan experiments are done for isoscalar targets, 
deuteron and carbon, and for tungsten which has a large neutron excess. The proton 
beam energy is 800 GeV. Forward production of the dileptons is dominated by beam- 
quark annihilation with a target antiquark in Eq. (|5.9| ). Therefore, the cross sections 
in the region xp > 0.1 can be used for investigating the antiquark distributions in 
the target. If sea-quark distributions in a nucleus are just the summation of proton 
and neutron contributions without nuclear modiflcation, the cross-section ratio for the 
tungsten (W) and the isoscalar targets (IS) becomes 



Ra{x) 



(7 1 six) 



1 + 



N — Z d{x) — u{x) 
A d{x) + u{x) 



(5.10) 



where A, Z, and A^ are the atomic weight, atomic number, and number of neutrons 
in the target nucleus. The dd annihilation is neglected because the uu dominates 
the cross section. Shadowing correction is applied to the tungsten data at x < 0.1 
in the following way. First, A-dependent shadowing factor ash is determined from 
EMC, NMC, and E665 shadowing data for D, C, and Ca. Then, the tungsten cross 
sections with x < 0.1 are corrected by a^ = a^A"'"'^. The obtained ratios are shown in 



Fig. Q^ together with three theoretical expectations. The solid curve is a pion-model 
prediction ||^, |B^ in section |4.3.1| . The dashed one is a simple parametrization in 
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section |2]^ for explaining the NMC data: d — u = A{1 — x)^ with A = 0.15(1 + b) 
and 6=9.6 [H^. The dot-dashed one is a chiral model result [^, ^ in section |4.3.2| . 
As it is obvious from the figure, the data do not reveal significant flavor asymmetry. 
Although the E772 data are consistent with the flavor asymmetric model predictions, 
they could be also explained by the flavor symmetric sea by considering experimental 
errors. The tungsten is a heavy nucleus, so that there could be a signiflcant nuclear 
effect on the flavor distribution (see section |4.6|) . The E772 collaboration also showed 
the xp distribution of p+d data in connection with the flavor asymmetry. However, 
there is also no evidence for the asymmetry. These results are somewhat in conflict 
with the NMC result and other Drell-Yan data. 

There are Drell-Yan data for various nuclear targets; however, the NA51 data point 



at a;=0.18 [^ is the only existing one which is extracted from the p-p and p-d Drell- 
Yan experiments. Considering the x range of the NA51 measurements, we neglect 
shadowing correction in the deuteron. Then, the asymmetry can be written as A^y = 
2a^^/a^'^ — 1. The NA51 collaboration used 450 GeV primary proton beam from the 
CERN-SPS. The targets are liquid hydrogen and deuterium. The accepted rapidity 
range is from —0.5 to 0.6, and the muon mass region M^^ >4.3 GeV is used for the 
analysis. The obtained asymmetry is 

Ady = -0 -09 ± -02 (stat.) ±0 .025 {sy St.) . (5.11) 

As the valence-quark value Xy = u^/d^ at x=0.18, they take Ay=2.2 averaged over 
the parton distributions, MRS-S'q, MRS-D'q, and MRS-D_', GRV-HO, and CTEQ-2M. 
From Eqs. ( |5.8|) and ( |5.11 ), the observed asymmetry becomes 



u 

-^ = 0. 51 ± 0. 04 ( Stat. ) ±0. 05 (sy St.) at x = 0.18 . (5.12) 

d 

This is a clear indication of the flavor asymmetry in the light antiquark distributions. 
There is an excess of d-quarks over u in the nucleon, and the NA51 result agrees with 
the tendency obtained by the NMC. Unfortunately, the only one data point at x = 0.18 
is available in the NA51. More complete Drell-Yan experiments at the Fermilab (E866) 
should give a clearer answer to the flavor symmetry problem [^. At this stage, the 
preliminary E866 data seem to show the NMC type asymmetry in addition. 

The E772 and NA51 data are compared with various model predictions in Refs. 



18, 43, 37, ^ 39, 70, 71, 72, 73 1. The theoretical works are done mainly to compare 
the mesonic calculations with the Drell-Yan data. At this stage, the mesonic models 
could be consistent not only with the NMC result but also with the E772 and NA51 
Drell-Yan data. 
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5.1.1 Fermilab-E866 results 

The Ferniilab-E866 collaboration reported high statistical experimental results for the 
ratio of the Drell-Yan cross sections crp(i/2o"pp |J20| . The dimuons are measured in 800 
GeV proton scattering on the liquid hydrogen or liquid deuterium target. Six month 
data were collected until March of 1997, and 350,000 Drell-Yan events were obtained. 
The goal of the experiment is to measure the cross-section ratios with 1% accuracy in 
the X range 0.03 < a; < 0.15. The accuracy becomes worse in the larger x region. At 
large x^, the ratio is approximated as [cTpd/2app]x >o ~ [1 + d{x)/u{x)]/2, so that it 
is possible to extract the distribution ratio d{x)/u{x) from the large x^ data of the 
cross sections. They compared the obtained ratios apd/la^p with the flavor symmetric 
distributions [CTEQ4M(m = d) and MRS(SO)] and asymmetric ones [CTEQ4M(m ^ J), 
GRV, and MRS(G)]. The symmetric CTEQ4M(m = d) curve is obtained by modifying 
the distributions a.s u = d = {u + d)/2. The accurate part of the data in the region 
0.03 < x < 0.15 tends to agree with the asymmetric distributions, in particular with 
the CTEQ4M('u ^ d). However, it is also interesting to find that the data seem to 
deviate from the present asymmetric parametrizations in the larger x region. Although 
their results are still preliminary, the accurate E866 data confirm the flavor asymmetry 
conclusions of the NMC and NA51. 
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5.2 W and Z production 

Instead of the virtual photon production in the Drell-Yan case, weak boson production 



could also have information on the antiquark distributions fT^, |78|, Q. There are 
existing CDF data for the charged lepton asymmetry (or W charge asymmetry) in the 
p + p reaction: pp — >• W^X -^ {i^i'()X [^. They are first analyzed in Ref. |^ 
in connection with the u — d distribution. Although the CDF data constrain the u/d 
ratio in the region of x = M^^ j ^ = 0.045, they are consistent with the symmetric sea 
u = d. However, the p+p reaction is not very sensitive to the sea-quark distributions as 
we discuss in this section. Therefore, future p + p colliders such as Relativistic Heavy 
Ion Collider (RHIC), rather than the p + p-, are crucial for investigating the flavor 
asymmetry in W and Z production processes. We discuss the sensitivity of W^ and Z^ 
production cross sections on the u/d asymmetry based on Ref. [|79| in the following. 

We show how the W production processes in the p + p collider could be used for 
probing the flavor asymmetry. The W^ production cross section is given by parton- 
subprocess ones together with parton distributions in the colliding hadrons p9|: 



a{p + p^W+X) = \ I dx, [ dx2y2q{xi,M^)q'{x2,M^)a{qq' ^W+) , 

^ ^0 Jo g,^, 

(5.13) 
where 1/3 is the color factor 3 • (1/3)^ = 1/3. The subprocess cross section is given by 

da{qq' -. W^) = (0^ ^ ^ l^^^^' ^ ^^^\' W S\Pi + P2 - p) ^^^ , 

(5.14) 

where pi, p2, and p are q', q, and W'^ momenta respectively, and s is given by s = 
(Pi +^2)^- The matrix element is 

Miqq'^W+) = -tV,,,-^ei\p)vip,)^ra-l5)uiP2) , (5.15) 

and the Cabibbo mixing is used in our calculation: Vud = cos 6'c, Vus = sin^c, Kd = 
— sin 6'c, and Vcs = cos 6c- Taking the spin summation, we obtain 

J2\M{qq' ^W-')\' = ^GfM^\V,,,\' , (5.16) 

pol 

with the Fermi coupling constant Gp/V^ = 5'^/(8M^). Noting 5^{pi + P2 — 
p)d^p/{2Ep) = 6{s — M^) and dxidx2 = dsdxp/[{xi + X2)s], we have the iy+ produc- 
tion cross section in the p + p reaction in terms of the parton distributions: 
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dxp o \xi -\- X2 / 

+ sin^ 6'c [u{xi)s{x2) + s(xi)m(x2)] } . (5.17) 

The dominant processes of producing W^ are u{xi)+d{x2) — > IV"^ and u{x2)+d{xi) — > 
iy+; however, the first one becomes much larger than the second at large xp- Therefore, 
the cross section is sensitive to the d distribution at large xp- On the other hand, the 
cross section for the W~ production is given in the same way: 



dap+p^w- _ v27r 



Gp [ ) { cos^ 9c [u{xi)d{x2) + d{xi)u{x2)\ 

dxp 6 x'^i ^ ^2 J 



+ ?,iY? 6c[u{xi)s{x2) + s{xi)u{x2)\] . (5.18) 



At large xp, it is sensitive to the u distribution instead of the d in the W^ case. This 
difference makes it possible to find the difference u — d. Because the Cabbibo angle is 
small, the sin^ 6^ terms are neglected for simplicity in the following discussions. The 
W^ production ratio is then given by 



dapj^p^w+ /dxp u{xi)d{x2) + d{xi)u{x2) 
dap^p^w- Idxp u{xi)d{x2) + d{xi)u{x2) 






At large xp (large Xi), the antiquark distribution g(xi) is very small, so that the above 
equation becomes 

«„(..» 0)«^^ . (5,20) 

which is directly proportional to the ratio d/u. 

The situation is very different in the p + p reaction case. Replacing the parton 
distributions in Eq. ( p.l9| ) by ^(0:2) — > ^(3^2) and q{x2) — > ^'(3^2), we obtain the ratio 

P / _ u{xi)d{x2) + d{xi)u{x2) , u{xi) d{x2) , . 

''-''^''^^~ u{x,)d{x2)+d{xMx2) "-»° d{xi)u{x2) ' ^ ' 

In the p + p reaction, the ratio is no more sensitive to the u/d asymmetry. How about 
the Xp ^ region? We find from Eq. ( ^.211) that the p + p ratio is independent: 



Rp+p{xp = 0) = 1, even though the p + p ratio is still sensitive to the flavor asymmetry 
at Xp = 0: Rp+p{xp = 0) = [u{x)/d{x)][d{x)/u{x)]. From these discussions, it is more 
appropriate to use a. p + p collider in finding the u/d asymmetry from W production 



data. This fact is numerically shown in Fig. pl5. 
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The ratios in the p + p and p + p reac- 
tions are evaluated at y/s=500 GeV in Fig. 
|5.5| by using various parametrizations for the 
parton distributions |]79[. The distributions 



^ife-- 


__EHQ 


p + p collisions ,^ 


_^> 


S^.,CTEQ2^ 


__^— -^^ 




■wz^^ _/ ' 






- ■ •^T^-^^MRSD-'' -- 


'a) 




--"""\\D01.1 ' 



p + p collisions 




are evolved to the scale Q^ = M^. The 
figures a) and b) show the p + p and p + p 
results respectively. The dashed curve indi- 
cates the results of using the flavor symmet- 
ric {u = d) DOl.l distributions. The others 
are the results for flavor asymmetric distribu- 
tions (MRSD_', CTEQ2pM, ES, EHQ). Be- 
cause the NA51 result ruled out the MRSD_' 
distribution, the small difference between the 
flavor asymmetric MRSD_' and the symmet- ^' 

ric DOl.l should not be taken seriously. As Figure 5.5: H/='= production ratios a) in p+p 
we expected, the p + p reaction is sensitive to ^^"^ ^) mp + p (taken from Rcf. |79[| ). 
the parton-distribution models, in particular 

the light antiquark flavor asymmetry, not only in the large \xf\ region but also in the 
Xi? ~ region. On the other hand, the p + p reaction is almost insensitive to the 
asymmetry. The model dependence appears only in the very small xp- 

The W production processes in the p + p and p + d reactions could also be used for 
studying the flavor asymmetry. The cross-section ratio is 

^ dapj^p^w+ /dxp u{xi)d{x2) + d{xi)u{x2) 



R'{xf) 



da„ 



(5.22) 



'p+d^w+ldxp u{xi) [u{xi) + d{x2)\ + d{xi) [u{x2) + d{x2)] 

by neglecting nuclear corrections in the deuteron. Although it is independent of the 
sea distributions at small xp'- R'{xf <C 0) = 1 -|- [u{x2) — d{x2)]/[u{x2) + d{x2)], large 
xp data are useful: 

d{x2) 



R\xf) 



U[X2J 



at Xf ^ 



(5.23) 



U{X2) + d{x2) 

In the similar way, Z^ production data in the p + p and p + d reactions are valuable. 
The Z^ production cross section in the p + p is 



da. 



dxp 



n 



3^2 



Gf 



X1X2 

Xi + X2 



^Xw + yxD [uixi)u{x2) + u{xi)u{x2)] 



+ (1 - i^Xw + -xi) [d{xi)d{x2) + d{xi)d{x2) + s{xi)s{x2) + s{xi)s{x2)] 



(5.24) 



where Xw is given by the Weinberg angle as Xw = sin^ 9w- Taking into account the 
dominant uu contribution at large xf, we obtain the ratio 



R"{xf) ^ 2 d^P+P^z^/d^^ 



da. 



p+d^Z 



o/dxF 



1 + 



U{X2) - d{x2) 
U{X2) + d{x2) 



a.t Xf ^ ^ 



(5.25) 
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We find that not only the W^ production but also the Z^ production could be 
used in determining the u — d distribution. Because p + p and p + d reactions are very 
sensitive to the asymmetry, future colliders such as RHIC should be able to find the 
antiquark asymmetry by the W^ and Z° production measurements. 



5.3 Quarkonium production at large Xp 

We discuss the possibility of finding the u — d distribution in quarkonium production 
processes. J/xp production data have been used in extracting gluon distributions in the 
nucleon and in nuclei. Because the dominant process is the gluon fusion gg ^ cc -^ 
J/ip, qq annihilation is in general a small effect. However, the qq process could become 
important at large \xf\. 

The mechanism of producing the quarkonium is a strong interaction, which makes 
the description more model-dependent than the electromagnetic Drell-Yan case. A 
popular description is a color-singlet (and recent color-octet) model, which includes 
99: 91: 91: and qq fusion up to a^. Instead of stepping into the detailed production 
mechanism, we discuss general features by selecting a simpler one, the semi-local duality 
model. The quarkonium production processes are analyzed in this model, and the 
results are related to the flavor asymmetry in Ref. [^ . 

The cross section for a QQ pair production is given by parton subprocess cross 
sections multiplied by the corresponding parton distributions 



dxpdr y/xp"^ + 4:t'^ 



G{xi)G{x2) a{gg —> QQ; m^ 



+ Yl {l'ixi)<l\x2) + q'{xi)q'{x2)} cT{qq ^ QQ]!!!"^ 

i=u,d,s 



. (5.26) 



The only gg and qq type subprocesses are taken into account in the above expression, 
and (T{gg -^ QQ; w?) and a{qq — >• QQ; rn?) are the corresponding cross sections. The 
variables xi and X2 are fractional momenta carried by the projectile parton and by the 
target one, and xp and r are given by xp = xi — X2 and r = mj y/s with the invariant 
mass of the QQ pair m. The subprocess cross sections are 



(y{.99 -" QQ;'m^ 



27m6 
7ra|_ 



(m^ + 2mi,) X 



{m^ + Am^rriQ + ml) ln( 



i (7m2 + 31ml) A 



m^ + A, 
m'^ — A' 



(5.27) 



where tuq is a quark mass and A is given by A = y w^'^ — Arn^mQ^. According to 
the semi-local duality model, the quarkonium production cross section is obtained by 
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integrating the subprocesses cross section from the QQ threshold to the open charm 
(beauty) threshold: 



da, 



dxp 



hrriQ/^ dxpdr 



(5.2J 



where F is the probability of a ^/ip (T) creation from the QQ state. 

We hope to find the antiquark flavor asymmetry from these quarkonium produc- 
tion processes. The gluon-gluon fusion process dominates the cross section in general. 
However, the qq processes could become more important in certain kinematical regions. 
The qq fusion contributions in the p + p collision are 

dap^p oc u{xi)u{x2) + u{xi)u{x2) + d{xi)d{x2) + d{xi)d{x2) ■ (5.29) 

In order to find the u — d distribution, the p + d reaction has to be studied in addition: 



dapj^d oc [u{xi) + d{xi)][u{x2) + d{x2)] + [u{xi) + d{xi)\[u{x2) + ^(^2)] 
where the isospin symmetry is assumed. The cross-section ratio 

da{p + p — > J /ip{T))/dxF 



R{xf) 



da{p + d —^ J/ip{T))/dxF 



(5.30) 



(5.31) 



should be sensitive to u — d particularly at large xp- From Eqs. ( |5.29| ) and (|5.30| ), it 
is obvious that the ratio is R{xf) = 1 in the flavor symmetric case u = d. Therefore, 
the deviation from unity is a signature of a finite u — d distribution. 

The ratio R{xf) is evaluated for the 800 
GeV proton beam in Fig. |5.(j| , where the up- 
per (lower) figure shows the J/ip (T) produc- 
tion results. Input parton distributions are 
the same in Fig. p.5| . The ratio is unity if the 
sea is flavor symmetric, and it is shown by the 
solid lines DOl.l in Fig. ^.6| . As we expected, 
effects of the flavor asymmetry become con- 
spicuous at large \xf\- This is because the 
gluon distribution G{x) is much smaller than 
the quark one q{x) at large x, and the cross 
section is dominated by the qq fusion pro- 
cesses. The parton-distribution dependence 
is more evident in the T production. Be- 
cause T is more massive than J/ip, we have 
(a;iX2)T ~ '^'rnB/y/s > 2mcl \fs. The T process is sensitive to the larger x region, so 



J/^ Production at 800 GeV 



r Production at 800 GeV 



Figure 5.6: J/V' and T production ratios 
(taken from Ref. 80 ) 



that the flavor-asymmetry effects become more conspicuous. 
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The above results show interesting contributions from the u/d asymmetry. There- 
fore, measurements of the J /ip and T production cross sections for the proton and 
deuteron targets at large xp should also be able to clarify the u — d distribution prob- 
lem. However, the semi-local duality model is probably too simple to explain the 
quarkonium-production processes. It is now revealed that the prediction of a more so- 
phisticated color-singlet model is inconsistent with the J /ip and ip' production data by 
the CDF. The discrepancy could be understood by the color-octet mechanism. There- 
fore, a better model analysis is necessary in comparing the theoretical results with 
future experimental data. 

5.4 Charged hadron production 

Semi-inclusive reactions in the electron or muon scattering could be used for finding 
the antiquark distributions. In particular, charged-hadron production could have in- 
formation on the flavor asymmetry. The hadron-production cross section is written 
by the lepton and hadron tensors in the same way with the inclusive one in section 
|2.1| . However, in spite of the fact that the semi-inclusive process is dominated by the 
light-cone region, the operator product expansion cannot be applied ||9^ . It is because 
the summation on X in the final state | X; p^, Sh > cannot be taken independently from 
the hadron state | ph, Sh >■ Therefore, we discuss the theoretical analysis on a relation 
between the charged-hadron-production cross section and the distribution u — d hj 
using a quark-parton model [|7^]. 

In the parton picture, the semi-inclusive cross section is given by [|[| 

1 da'^^{x,z) _Z.eU.{x)Dl{z) ^ ^^^^^^ 



(^n{x) dz Y^i^^lfii^) 

where fi{x) is the quark distribution with flavor i and momentum fraction x, and D^{z) 
is the z-quark to /i-hadron fragmentation function with z = Eh/v. The numerator for 
charged hadron production is 



N^'^^Y.'U\{^)D 



^^'z) 



= l^D^ + \^Dt + \dD^+^-dD^^ + ^-sDf + YsDf . (5.33) 

Assuming the isospin symmetry in the parton distributions, we consider a combination 
of proton and neutron cross sections: 



R{x, z) 



{NP+ - Ar"+) + {NP- - A^"-) 
{NP+ - A^"+) - {NP- - iV"-) 
u{x) - d{x) + u{x) - d{x) 4D+(z) +ADt{z) - Dj{z) - D^{z) 



u{x) - d{x) - u{x) + d{x) 4:D+{z) - 4D^(z) - DJ{z) + DJiz) 
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. (5.34) 



Here, A^^"*" {Np~) and A^""*" ( A^"~) correspond to the production processes of positively 
(negatively) charged hadrons from the proton and the neutron respectively. If the 
denominator and numerator are integrated over x individually, the Gottfried sum is 
obtained from the numerator integral, and the denominator becomes a sum for the 
valence quarks: 

Jdx{ {NP+ - A^"+) + {NP' - A^"-) } 



J dx { {NP+ - Ar"+) - {NP- - N^ 
^^ ADUz)+^Df{z)-Dj{z)- 


-)} 


^'^4D„n.)- 


-4Df{z)-DJiz) + D±{z) 



(5.35) 

According to this equation, if the Gottfried sum rule is violated, it should appear 
in the charged-hadron-production asymmetry. Available EMC data WM are analyzed 



by using Eq. (|5.35|) [[7^] . Contributions from pion, kaon, and (anti)proton production 



processes are taken into account in evaluating the fragmentation function, for example 

D: = D:^+D^^ + DP . (5.36) 

Isospin and charge conjugation invariance reduces the number of fragmentation func- 
tions for the pion: 

D^Df=Df = Dl=Df , 

D^Df =Df = D:- =Df . (5.37) 

In the kaon case, the reflection symmetry along the V-spin axis {D^ = D^ ) is used 
in addition to the isospin and charge conjugation invariance: D^"" = D^ = D^ , 
f)K ^ ^K~ ^ ^K+^ f)'K ^ jjK+ ^ jjK- ^ JJK+ ^ jjK' Furthermore, D^ and D'^ 

are assumed equal. Similar equations are taken for proton and antiproton production: 
DP = DP^ = Dl = Dl = D? DP = Dl = Dl = Dp = D^. Experimental information 
is provided for these fragmentation functions. In particular, we use parametrizations 



fitted to the EMC data 0: 
b{z) 1 



D{z) l + z ' 

D^iz) b^iz) 1 - z 

— 1^ = 0.35^ + 0.15, — ^ = QAhz , 

D{z) ' D{z) l + z ' 

DPiz) DPiz) l-z , , 

7J^ = 0.20, _l).o.l2— . (5.38) 

With these experimental parametrizations for the fragmentation functions, the ratio of 
the cross sections Q{z) becomes 
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Experimental data are given for 
{da^/dz) /aN, so that Fi{x) is multiplied in 
getting A^^'^. The experimental data zQ'^^ 
obtained in this way are compared with Eq. 
( |5.39| ) in Fig. |5]^. The upper curve is ob- 
tained by assuming the Gottfried sum Iq = 
1/3 in Eq. (|5.39| ). On the other hand, the 
hatched area is based on the 1991 NMC re- 
sult. It is interesting to find the difference 
between the two results in the semi-inclusive 
processes. However, as it is obvious from the 
figure, we cannot judge whether or not the 
sea is u/d symmetric from the data. 

The recent HERMES preliminary data seem to be accurate enough to find the u/d 

and vr" 



Figure 5.7: Charged- hadron- production ratio 



zQ'''^{z) (taken from Ref. ]76|) 



asymmetry 

-R(x, z) for the pion by 



17|1 . The following vr 



production ratio is related to the function 



r(a;, z) 



JSJP-K- 



AT" 



Rn(x, Z) — 1 



(5.40) 



The obtained data of r(x, z) in the range 0.1 < a; < 0.3 agree well with the NMC flavor 
asymmetry, and they are significantly different from the symmetric expectation. The 
HERMES results will be submitted for publication in the near future. 



5.5 Neutrino scattering 

Neutrino interactions are useful for determining the valence-quark distributions by tak- 
ing advantage of parity-violation terms. On the other hand, neutrino-induced dimuon 
data are used for determining the s and s distributions, so that the neutrino interac- 
tions could be valuable also for determining the light antiquark distributions u and d. 
We discuss what kind of cross-section combination is appropriate for finding the u — d 
distribution. The neutrino reaction via the charged current is given by the amplitude 

is 



Gf/V2 



1 + Q-'/M^ 
so that the differential cross section becomes 



da 



2 u{k')r{l - ^,)u{k) < X I j;-^(0) I p, a > , (5.41) 



M Gl 



s - M2 (27r)2 (1 + Q^/Ml,Y '^^ E' 



The leptonic tensor is given by 

^'^' = E,.v [^(fc',A')7^(l-75)n(A;,A)r [u{k\>!)Y{l - ^,)u{kA)\ 
= 2 ( k^'k"' + k'^k" - k ■ k'g^"" + ie^'^'P^kX ) , 



(5.42) 



(5.43) 
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where £^^^p'' is an antisymnietric tensor with e^^'^^ = +1. The last term does not appear 
in the electron or muon scattering because it is associated with the parity violation in 
weak interactions. This term makes it possible to probe new structure in the target 
hadron. There exists an antisymmetric term under the /x ^-> i/ exchange in addition to 
the hadron tensor in Eq. (127 



- ljW,F_^,,^jfq' . (5.44) 

The W'i structure function is proportional to the difference between left- and right- 
transverse cross sections for the W boson. With these structure functions, the cross 
section becomes 



da^ GlE'^ 



dQdE' 27r2 (1 + g7M2 ^2 



'WJ 



,^, o & E -\- E ^^, o 



2 



2 Wi sin^ - + W2 cos^ - T —j^ W3 sin^ - 



(5.45) 



where ± indicates W^ in the reaction. Structure functions Fi, F2, and F3 are defined 
by Fi = MWi, F2 = UW2, and F3 = UW3. 

On the other hand, the charged-current process is described by neutrino-quark 
interactions with the current 

Jf, = u{x) 7^ (1 - 75) [ d{x) cos 6c + s{x) sin 6^ ] 

+ c{x) 7^ (1 — 75) [ s{x) cos 9c — d{x) sin 9c ] ■ (5.46) 

Comparing a calculated cross section in the parton model with Eq. ( |5.45| ), we express 
the structure functions in terms of quark distributions 

Fi = F2/2x , 

F^^ = 2x{d+s + u + c) , 

F^'P = 2x{u + c + d + s) , 

xF^^ = 2x{d+s-~u — c) , 

xFl'P = 2x{u + c-d-s) . (5.47) 

For the time being, we discuss only LO contributions without NLO corrections from the 
coefficient functions. Because the antiquarks have negative parity, there are negative 
signs in the F3 structure functions. Combining the v and v F^ structure functions, we 
obtain the valence quark distribution [F^^ -|- F^^)/2 = u^ + d^ with the assumptions 
s = s and c = c. Therefore, it is the advantage of neutrino reactions that the valence 
distribution can be determined. However, they are also used for studying antiquark 
distributions. In fact, neutrino-induced dimuon data enable us to determine the s 
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From Eq. ( p.47| ), it is possible to combine 
the F2 and F3 structure functions for the pro- 
ton and the deuteron in order to get the flavor 
asymmetry 



31 



u-d=- {Fl^'lx - F, 



3 ) 



1 



- {F^^/x - F. 



vd\ 



(5.48) 



Figure 5.8: Cont ribut ions to t he 



by neglecting nuclear effects in the deuteron. It 

may also be defined in terms of Fi and F3 struc- 

ture functions 22 
distribution in Eq. ( 5.48 ) or ( 5.49 ) (a) from ^—^ 

quarks and (b) from antiquarks. 1 _ 1 

u - d = -{F^ - F^) - -{F^ - Fl^) , 

(5.49) 

if the s and c distributions can be neglected! 
These expressions are correct only in the LO. However, they are practically convenient 
because the charged reactions are easier to be measured experimentally. Unfortunately, 
present data are not accurate enough to be used for finding the u — d distribution. To 
be precise, Eqs. (|5.48|) and ( |5.49|) are not appropriate in defining the u — d distribu- 
tion if NLO effects are taken into account. For example, the coefficient functions are 



different in Fi and F3 structure functions, so that [F^^ — F^^^jl — [F^^ 



Fl^)l^ has 

a contribution, already of the order of a^, from the quark in Fig. p.8| (a) in addition 
to the antiquark in Fig. |5.8| (b). Although it is useful in getting experimental informa- 
tion o\iu — d from Fx and F3 in the charged reactions, it is not a precise definition if 
higher-order corrections are taken into account. 

A consistent way is to use both the charged and neutral current reactions |2^ . The 
neutrino-quark interaction via the neutral current is described by [Q 



'^'' = X] 9 ^*(^) ^^ {9Li{^ - 75) + ^ra(l + 75) ] (liiA 



(5.50) 



where gu and gB.% are defined by the Weinberg angle Qw 

4 . .„ 4 

3 



9Li 



1 - - sin^ e 



w 



dm 



2 
3 



- sin^ Ow 
2 



— l + -sin^6'iy , = -F- sin^6'vi/ 



for i=u, c 



for i=d, s 



(5.51) 



With the neutral-current observables, the distribution u ~ d can be defined by only 
one type of the structure functions, for example Fi. The neutral-current structure 
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functions Fi for neutrino and electron reactions become 

1 2 



F^^ = g sec Ow 



4 



- sin Ow + Q sin 9w 



1 2 

- sin^ 6'vi/ H — sin^ 6^ 



{m(x) + u{x)} 
{d{x) +d{x)} 



w 



n 



ep^eX 



g"^ sin^ 9w 



4 1 - 

- {u{x) + u{x)} H — {(i(x) + d{x) } 
9 9 



(5.52) 



(5.53) 



According to the definition of Ref. [Q, it is given by combining the charged and 
neutral current Fi structure functions as 

Fi.^) = 2 [-^i(^) ~ {Fi{,x)^p^e+x - Fi{x)yp^e-x}] = u{.x) - d{x) , (5.54) 

where Fi{x) is defined by the structure functions in the neutral current reactions: 



Fi(x) 



(I sin'^e^. - I) secW^ 



-FAx] 



up~*uX 



9 (I - stn^e^ + f stn^e^) 
(I sin'^e^ - I) sin'^e^ 



Fi{x)ep^eX ■ 

(5.55) 



Because the s and c distributions are neglected in the above discussion [^ and in Eq. 
(|5.49|) , it is necessary to subtract out these contributions by combining Eqs. (|5.54| ) and 
( |5.49| ) with the deuteron Fi structure functions. 

It is impossible to obtain the u/d asymmetry from present neutrino data. How- 
ever, we hope that much better data will enable us to extract the flavor asymmetry 
distribution. 



5.6 Experiments to find isospin symmetry violation 



We discussed in section |]5| that the violation of the Gottfried sum rule could be due to 
the isospin-symmetry violation instead of the flavor asymmetry. These two mechanisms 
cannot be distinguished at this stage. The isospin-violation effects are believed to be 
very small in the structure functions. However, it is important to confirm this common 
sense experimentally. Various processes are discussed in the following for finding the 



effects of the isospin-symmetry violation in the antiquark distributions [50 



The F2 structure functions in neutrino interactions are useful in distinguishing 
between the two mechanisms [5^] . The difference between proton and neutron structure 
functions is 

dx 1 



usv 



X 2 



K^(x) + K^(x)-Fr(a:)-K 



dx [ { u{x) + d{x) + s{x) + c{x) }p — { u{x) + d{x) + s{x) + c(x) }„ ] 

(5.56) 
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If the failure of the Gottfried sum is entirely due to the flavor asymmetry, the integral 
vanishes Ijsv = 0- On the other hand, if it is entirely due to the isospin violation and 
if the s and c terms can be neglected, the integral is Ijsv = —0.336 ± 0.058. Because 
the flavor asymmetry does not contribute, the sum Ijsv should give a clue in finding 
an isospin-violation sign. 

Isospin- violation effects on the Drell-Yan processes are also discussed in Ref. [ ^ . 
In the pion scattering case n^A —>■ i~^i~X, we consider the difference of nuclear cross 
sections at large Xj,: 

_ 4 [aJTr+A,] - aJTr+Ao) ] + [aJTr-A,) - a{n~Ao) ] 

The Aq and Ai denote different nuclear species, but we may choose Aq as an isoscalar 
nucleus and Ai as a neutron-excess nucleus. With the isospin symmetry assumption, 
it becomes 

fi„-= '°'^' ";'<""'' . (5.58) 

where e is a neutron excess parameter e = N/A — 1/2. On the other hand, if the sea 
is flavor symmetric with the isospin violation, the Drell-Yan ratio becomes 

_ 50(e,-eo)(g^-g") 

where the isospin symmetry is assumed for the valence-quark distributions. Similar 
equations are obtained for proton Drell-Yan cross sections. The p-n cross section 
asymmetry is given in the isospin symmetry case as 



Ady 

(5.60) 



{Auy — dv){u — d) + {uv — (iy)(4-u — d) 



{Auv + dv){u + d) + {uv + dv){Au + d) 

On the other hand, it is given in the isospin-violation case as 

^ _ {Auy - dy) 5 {qP - g")/3 + {uy - dv){qP + 8g")/3 
^^ 9 (o-PP + crP'^) ■ 

The details of the Drell-Yan cross sections and the asymmetry are discussed in section 



0| . From these equations, we find that the Drell-Yan cross sections could be interpreted 



in principle either by the flavor asymmetry or by the isospin violation. Both effects 



are taken into account to explain the NA51 result in Ref. [^. The obtained result 
indicates that the ratio u/d could be larger than the NA51 value at the cost of isospin 
symmetry violation. However, it is not possible to separate these two contributions 
clearly. 
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It is shown in Ref. [^ that the flavor asymmetry could be found in semi- 



inclusive leptoproduction of charged hadrons. The number of produced h hadrons 
in the lepton-nucleon scattering at Bjorken x and z = Eh/v is given by N^^{x,z) = 
Si e^ qf {x) D'^z), where Df is the fragmentation function. The details of the charged- 



hadron production are discussed in section |5^ . The following equation is obtained for 
flnding the flavor asymmetry: 



Q{z) 



3 J 0-50^^ + 3.1Z + 7.6 

^ 3.2^2 + 112 + 0.84 ' ^ ' ' 



where the notations A^^+, A^"+, A^^ , and A^" are the same as those in Eq. ( |5.34| ). 



The isospin symmetry is assumed in the above equation. If the sea is flavor symmetric 
and if the isospin symmetry is violated, the above quantity becomes [pOl 



_ 4[Dt{z) + DUz)]il-26q) - [Djjz) + D±{z)]il + 26q) 
^^'^ ^DUz)-Dt{z)]- [Djiz)-D±iz)] 

^3 0.80.2 + 3.37. + 7.63 

3.2.2 + 11^ + 0.84 ' ^ ^ 

where 6q = q^ — q"-. Both expressions have different . dependence, so that we should be 
able to distinguish the mechanisms if experimental data are accurate. At the present 
stage, charged-hadron production data are not accurate enough for flnding the discrep- 
ancy. 

The experimental studies would be difficult because the isospin effects are consid- 
ered to be very small theoretically. However, accurate experimental data are desper- 
ately needed in order to shed light on the isospin-symmetry violation in the antiquark 
distributions. 
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6 Related topics on antiquark distributions 

As a topic of flavor asymmetry, we have discussed the hght-antiquark-distribution dif- 
ference in the nucleon. There are other important issues on the antiquark distributions. 
We briefly comment on related topics. 

First, the u/d asymmetry in hyperons could be studied if charged hyperon beam 
becomes available in future. For example, a possibility to find the asymmetry in S^ 
is investigated in Ref. |100]. In a naive quark model, they consist of S^(mms) and 



Yi^[dds). The Pauli-blo eking and meson-cloud models predict d excess over u in S+ 
and u excess over d in S~. It is an interesting test of the theoretical models in section 
0. The Drell-Yan cross section for the S+p reaction at y = is given by 






9^7 



4 _ _ 1 - 

- { Up{x)u-E{x) + UT.{x)up{x) } + t: { Up{x)dY.{x) + SY.{x)sp{x) } 

y y 



(6.1) 



where only valence-sea annihilation terms are retained, and gs denotes the distribution 
in S+ (gs = i?E+)- In calculating cross sections for other reactions S^n, S+n, and S~p, 
we assume isospin symmetry, Up = dn, Up = dn, uy.+ = d^-, uy.+ = (is-, together with 
the assumption ss+ = ss-- From the Drell-Yan cross sections with S^ beams on the 
proton and deuteron targets, we take the ratio 



R{x) = 






_ Mx) [rpjx) - rp{x)] - [1 - rp{x)rp{x)] , ^ 

" 5 {Tp{x) - 1] ^*^-" ' ^^-^^ 

where r^ = Up/dp, fp = Up/dp, and fs = us/^s- In this way, if r^ and fp are known 
from other experiments, fg could be measured by the hyperon Drell-Yan experiments. 
Second, we mentioned the s-quark distribution difference from the {u + d)/2. It has 
been measured experimentally by the neutrino induced opposite-sign dimuon events. 
This topic is studied within the meson-cloud models. For example, because the pions 
do not contain the valence s quark, their contributions to s and {u + d)/2 in the proton 
are different. It is particularly important in discussing the size of the nNN form 
factor and its relation to nuclear potentials. This topic is also discussed in the chiral 
field theory, so that the interested reader may look at the meson-model papers in the 
reference section. 



Third, the difference between s and s is also important ||101|| . Because there is no net 
strangeness in the proton, the integral of the difference has to vanish: J dx{s — s) = 0. 
However, x dependence of both distributions could be different. In fact, the proton 
virtually decays into for example K~^{us)A{uds), K~^ {us)T,^ {uds) , and K*^ {ds)T,^ (uus) . 
Within the three decay modes, the valence s is contained in the kaons and s is in 
the hyperons. Because the hyperon masses are larger than those of the kaons, the s 
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distribution is distributed in the outer side. It means that the s distribution is softer 
than that of the s-quark one. Of course, the s and s distributions should be dominated 
by the perturbative contributions. However, these could be canceled out by taking the 
difference s — s. It is impossible to find this kind of small effect at this stage 0. We 
hope to have much accurate data in future. 

Fourth, flavor asymmetry in polarized antiquark distributions should become an 
exciting topic in the near future. As far as the model is concerned, we have explained 



the flavor dependence in section [4.7| . However, because we do not have a variety of 
polarized data at this stage, it is very difficult to find the difference between it, d, and 
s distributions from experimental data. In any case, there is an attempt to study the 



fiavor decomposition by including semi- inclusive data in Ref. [ [L02|| . Future experi 



mental programs for the polarized fiavor asymmetry are for example the RHIC-SPIN 
and the Common Muon and Proton Apparatus for Structure and Spectroscopy 



(COMPASS) [|104|| . In the similar way with the unpolarized case in section |5.2| , the 
W^ production measurements by the RHIC-SPIN collaboration should enable us to 
find Au and Ad distributions. The strange polarization and other polarized valence 
and sea distributions will be measured in semi-inclusive reactions by the COMPASS 
collaboration. Much progress is expected on the fiavor dependence of the polarized 
antiquark distributions in the next several years. 

Fifth, there is a similar sum rule to the Gottfried in the spin-dependent structure 
function 61 for spin-one hadrons. This new structure function is related to quadrupole 
structure of the spin-one hadrons. Its sum rule was proposed in Ref. [ |lU5[j as 



/5 t 
dxbi{x) = Jim -3 ^]^^q(^) + ^Qsea , (6.3) 

where Fqlt = 0) is the quadrupole moment in the unit of e/M^ for a spin-one hadron 
with the mass M. The second term 6Qsea is the sea-quark tensor polarization defined, 
for example, 5Q£„ = / dx[86u{x) + 26d{x) + 6s{x) + 6s{x)]^ /9 for the deuteron. The 
distribution 6q is given by Sq = [g*^ — (g"*"^ + g~^)/2]/2, where the superscript indicates 
the hadron helicity in an infinite momentum frame. The Gottfried sum 1/3 corre- 
sponds to the first term lim^^o ~3iipFQ{t) = 0. Because the valence-quark number 
depends on fiavor, the finite sum 1/3 is obtained in the Gottfried. However, it does 
not depend on spin, so that the first term vanishes in the 61 case. The second term 
in Eq. ( |6.3|) corresponds to J dx{u — d) in Eq. (|2.12| ). Therefore, a deviation from 
the sum J dxhi (x) = should suggest the sea-quark tensor polarization as the Got- 
tfried sum rule violation suggested the finite u — d distribution. Recent studies indicate 
that the diffractive-nuclear-shadowing and pion-excess mechanisms produce a tensor 



polarization, which leads to violation of the hi sum rule |106 
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7 Summary and outlook 

The light antiquark distributions u and d had been assumed equal for a long time. 
The Gottfried sum rule can be derived with this assumption. Even though there were 
some experimental efforts to test the sum rule and the flavor asymmetry u — d, it 
was not possible to draw a reliable conclusion. However, recent accurate experimental 
measurements made it possible to find the difference between the u and d distribu- 
tions. The NMC finding of the Gottfried-sum-rule violation and the u/d asymmetry 
motivated us to study theoretical mechanisms and different experimental possibilities. 
The flavor asymmetry is now confirmed by the NA51 Drell-Yan experiment, and it 
is also suggested by the preliminary HERMES and E866 data. On the other hand, 
future experimental facilities should be able to pin down the u and d distributions. For 
example, Drell-Yan and W^-production measurements at RHIC should be very useful. 
In testing the Gottfried sum itself, we need to accelerate the deuteron at HERA. 

On the theoretical side, the perturbative corrections to the sum are very small. 
Therefore, the violation should be explained by a nonperturbative mechanism. Within 
the proposed models, the mesonic model is a strong candidate in the sense that it can 
explain the major part of the violation. The Pauli blocking effect is smaller than that 
of the mesonic model according to the naive counting estimate. Furthermore, if the 
antisymmetrization is considered in addition to the Pauli principle, both mechanisms 
could produce a u excess over d. Because there are other theoretical candidates as 
explained in this paper, we should investigate more details of these models in order to 
find a correct explanation. The flavor asymmetry studies provide us an important clue 
to understand nonperturbative aspects of nucleon substructure. Future experimental 
and theoretical efforts on this topic are important for understanding internal structure 
of hadrons. 
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